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INTRODUCTION  AND  SUMMARY 

1.1  PURPOSE  AND  OBJECTIVES  OF  THE  STUDY 

The  Navy-wide  component  of  the  current  DOD  Plan  for 
Atmospheric  Transmission  Research  and  Development  is  en¬ 
compassed  in  the  E-0  MET  Program  managed  by  the  Naval  Ocean 
Systems  Center  (NOSC) .  The  Navy's  broad  area  of  responsi¬ 
bility  within  E-0  MET  is  measurement  and  modeling  of  atmos¬ 
pheric  propagation  conditions  as  affected  by  the  marine 
environment.  Specific  E-0  MET  goals  include: 

•  Establishment  of  a  marine  aerosol  measurements  data 
base  including  vertical  distributions. 

•  Development  of  relationships  between  meteorological 
measurements  and  aerosol  models. 

•  Performance  of  meteorological  and  propagation  mea¬ 
surements  in  concert  with  systems  tests. 

•  Development  and  validation  of  systems  performance 
models  based  on  weather  climatology. 

A  key  element  of  the  E-0  MET  plan  is  an  E-0  systems 
performance  assessment  system  which  goes  under  the  name 
Prediction  of  Performance  Range  for  E-0  Systems  (PREOS) . 
When  implemented  on  a  desktop  Programmable  Calculator 
(PROCAL)  the  system  is  known  as  PROCAL-PREOS  [1] .  The 
primary  function  of  PROCAL-PREOS  is  to  provide  real  time 
shipboard  assessment  of  airborn  FLIR  performance.  As  a 
requirement  of  the  development  process,  PROCAL-PREOS  must 
be  designed  and  validated  against  the  expected  range  of 
the  operational  marine  environments.  It  is  to  this  broad 
objective  that  the  work  described  in  this  report  responds. 
Herein,  the  development  and  implementation  of  a  LIDAR  spe¬ 
cifically  for  the  E-0  MET  requirements  is  evaluated. 
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This  evaluation  has  taken  the  form  of  a  study  in  which 
several  distinct  but  nevertheless  key  areas  pertaining  to 
thj  development  and  deployment  of  a  shipboard  LIDAR  system 
have  been  analyzed.  It  is  intended  that  the  use  of  such 
a  LIDAR  will  provide  a  real-time  or  near  real-time  measure¬ 
ment  of  infrared  aerosol  extinction  coefficient  profiles, 
which  information  can  in  turn  be  used  as  input  into  a  real¬ 
time  FLIR  performance  prediction  such  as  PROCAL-PREOS  oper¬ 
ating  aboard  ship.  Such  a  LIDAR  measurement  system  is  in¬ 
tended  to  replace  the  use  of  meteorological  visual  range 
(visibility)  measurements  performed  by  human  observers, 
coupled  with  models  for  wavelength  scaling  of  these  data, 
and/or  the  sole  reliance  upon  predictive  models  which  use 
measured  wind  and  relative  humidity  values  to  estimate  in¬ 
frared  aerosol  extinction. 

Included  among  the  major  areas  of  investigation  in  this 
work  are  a  review  of  marine  aerosol  modeling  techiques  and 
comparisons  to  recently  measured  optical  and  aerosol  counter 
data  collected  in  marine  environments.  The  objective  of 
this  portion  of  the  study  is  a  realistic  characterization 
of  the  marine  environment  in  which  a  shipboard  LIDAR  will 
operate  so  that  meaningful  estimates  of  performance  can  be 
made  and  that  system  design  trade-offs  can  be  analyzed. 

These  include  selection  of  operating  wavelengths,  laser 
source  operating  mode  and  cw  output  power  or  pulse  energy 
and  transmitter  and  receiver  optical  system  size. 

Various  approaches  to  the  use  of  LIDAR  systems  for 
the  measurement  of  atmospheric  aerosol  attentuation  have 
been  proposed  and  tested.  Several  approaches  using  cw 
laser  sources  which  somewhat  relax  the  eye  safety  considera¬ 
tions  associated  with  the  use  of  high-peak-power  pulsed 
lasers  were  reviewed  in  this  study.  Some  of  the  more  im¬ 
portant  considerations  concerned  with  the  extraction  of 
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atmospheric  extinction  coefficients  from  pulsed  LIDAR  re¬ 
turn  data  have  been  considered  and  the  implications  for  use 
of  data  collected  at  shorter  wavelengths  for  the  estima¬ 
tion  of  aerosol  extinction  coefficients  at  longer  wavelengths 
have  been  considered  to  determine  the  optimum  choice  of 
system  operating  wavelength  and  the  inherent  limitations 
associated  with  this  choice. 

The  important  area  of  diverse  approaches  to  LIDAR  data 
signal  processing  was  also  reviewed  in  this  study  in  order 
to  best  satisfy  the  somewhat  conflicting  requirements  for 
rapid  data  reduction  and  the  generation  of  near-real-time 
information  versus  a  requirement  to  collect  and  efficiently 
process  and  manipulate  LIDAR  data  in  a  more  nearly  "raw- 
data"  condition,  so  that  important  but  lacking  information 
on  inversion  method  validity  and  suitability  can  be  further 
investigated. 

The  major  factors  involved  in  the  implementation  of 
a  LIDAR  measurement  system  in  a  shipboard  context  were 
considered  to  take  into  account  some  of  the  many  practi¬ 
cal  considerations  which  will  ultimately  have  a  major  im¬ 
pact  on  the  degree  of  success  with  which  such  a  system  can 
be  designed,  constructed  and  operated  reliably. 

Performance  estimates  for  an  exsisting,  prototype, 
hand-held  LIDAR  system  (the  US  Army  Atmospheric  Sciences 
Laboratory  Model-1  visioceilometer )  have  been  performed 
and  compared  to  ones  performed  for  an  upgraded  system 
which  would  use  a  more  energetic  laser  source  together 
with  a  larger  optical  system.  Performance  comparisons  of 
systems  operating  at  1.06  ym  to  ones  operating  in  the  3-5 
ym  and  8-12  ym  bands  were  performed  to  provide  guidance 
for  future  system  development. 
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1.2  STATEMENT  OF  TASKS 


The  several  wide-ranging  topics  included  in  the  pre¬ 
sent  study  were  divided  into  four  major  tasks  which  are: 

(1)  A  definition  of  the  marine  shipboard  environment, 

(2)  A  review  of  LIDAR  measurement  system  approaches 
and  data  processing  techniques, 

(3)  A  survey  of  exsisting  MIL-SPEC  and  ruggedized 
laser  and  E-0  hardware  which  might  be  incor¬ 
porated  into  a  shipboard  LIDAR  system  as  cri¬ 
tical,  key  building  blocks,  and 

(4)  Preparation  of  a  summary  of  findings  and  the 
development  of  recommendations  relating  to  the 
design,  construction  and  use  of  a  proposed 
shipboard  LIDAR  system. 

The  organization  and  content  of  this  report  closely  follow 
the  division  of  the  study  into  the  above  four  tasks. 

1.3  INTERIM  PROGRESS 

1.3.1  INTERIM  REPORT  #1 

An  interim  progress  report  [2]  on  the  present  study 
was  prepared  and  supplied  to  NOSC  in  September  1981.  The 
contents  of  that  report  have  been  substantially  incorpor¬ 
ated  into  this  final  report. 

1.3.2  OPTICAL  SOCIETY  OF  AMERICA  PRESENTATION 

A  verbal  presentation  describing  results  and  tenta¬ 
tive  conclusions  reached  in  the  study  and  available  at  the 
time  of  the  presentation  (October,  1981)  was  made  by  one 
of  the  authors  (JAD)  at  the  1981  Optical  Society  of  America 
National  Conference  held  in  Orlando,  Florida.  Copies  of 
the  viewgraphs  used  in  that  presentation  are  included  as 
an  Appendix  to  this  report. 
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1.4  ORGANIZATION  OF  THE  REPORT 


This  final  report  is  organized  into  five  major  sec¬ 
tions,  the  first  of  which  is  this  Introduction.  Features 
of  several  marine  aerosol  models  currently  in  use  and 
several  examples  of  recent  experimental,  infrared  aerosol 
extinction  data  are  discussed  and  compared  in  Section  2 
entitled  Characterization  of  the  Marine  Environment.  Mea¬ 
surement  approaches  and  techniques  for  both  pulsed  and  cw 
LIDAR  systems  are  presented  and  discussed  in  Section  3. 
Several  important  factors  involved  in  the  application  of 
a  LIDAR  measurement  system  to  a  shipboard  installation,  in¬ 
cluding  eye  safety  considerations  and  the  utilization  of 
MIL-SPEC  E-0  hardware  components  are  discussed  in  Section 
4.  Section  5  contains  a  statement  of  specific  recommenda¬ 
tions  which  have  resulted  from  the  review  and  analysis 
performed  during  this  study. 

1.5  SUMMARY 

A  study  has  been  performed  for  the  NOSC  to  review  and 
evaluate  several  factors  pertaining  to  the  Definition  of 
a  Shipboard  Visibility  Measurement  System.  This  study  was 
divided  into  four  major  tasks  which  include  a  characteriza¬ 
tion  of  the  marine  environment,  a  review  of  LIDAR  measure¬ 
ment  system  approaches  and  techniques,  a  survey  of  exsisting 
MIL-SPEC  and  ruggedized  E-0  hardware,  and  the  preparation 
of  summary  findings  and  the  generation  of  recommendations 
for  the  development  of  a  shipboard  LIDAR  system. 

The  major  conclusions  reached  in  this  study  are  as 
follows : 

1.  Regarding  Aerosol  Models 

•  The  models  predict  a  very  wide  range  of  visible 
and  infrared  aerosol  extinction  as  a  function 
of  windspeed  and  relative  humidity. 
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•  Infrared  aerosol  extinction  predictions  of  the 
Katz-Ruhnke  (K/R)  aerosol  model  are  lower  than 
LOWTRAN  5  for  5  m/s  windspeeds  and  higher  than 
LOWTRAN  5  for  windspeeds  >  15  m/s. 

•  Values  of  the  backscatter  coefficient  8  for 
1.06  urn  wavelengths  range  over  a  factor  of  300 
for  a  2  km  visibility  fog  compared  to  a  5  km 
visibility  oceanic  aerosol. 

•  The  ratio  of  8  to  the  aerosol  extinction  coef¬ 
ficient,  a  decreases  with  increasing  wavelength 
by  a  factor  of  30  for  wavelengths  from  4  ym  to 
10  ym. 

2 .  Regarding  Aerosol  Measurements 

•  Aerosol  data  show: 

-  no  clear  choice  of  model 

-  moderate  correlation  of  visible,  1.06  ym  and 
3-5  ym  scattering 

-  relative  independence  of  visible/near  IR  and 
10.6  ym  scattering 

•  Data  for  the  ratio  of  infrared  to  visible  aero¬ 
sol  extinction  coefficients,  a.,  /a  cc>  show 
large  amounts  of  scatter  and  usually  fall  be¬ 
low  predicted  values. 

3.  Regarding  LIDAR  System  Performance  Evaluations 

•  The  ASL  Model-1  visioceilometer  system  appears 
limited  to  operation  for  ranges  of  about  1-2  km. 

•  An  upgraded  1.06  ym  system  with  a  pulse  energy 
of  10  xJ/pulse  and  a  receiver  aperture  diameter 
of  30.5  cm  (12  inches)  promises  good  perfor¬ 
mance  for  ranges  >  5  km. 

•  Large  optics,  higher  pulse  energy  (>  10-1  J/ 
pulse)  and  heterodyne  detection  will  be  re¬ 
quired  for  10.6  ym  systems  due  to  decreased 
values  of  8  and  8/a  at  this  wavelength. 

•  The  condition  required  for  use  of  the  "Klett 
method"  for  LIDAR  data  inversion,  namely  that 
8  =  A  k  where  A  and  k  are  constants  and  not 
wavelength  dependent  is  not  borne  out  by  calcula¬ 
tions  of  the  ratio  of  scattering  to  extinction 
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at  10.6  ym,  rendering  the  use  of  this  approach 
at  10.6  ym  questionable. 

•  A  reliable  method  for  estimation  of  a.,,  from 
analysis  of  1.06  ym  LIDAR  returns  has  not  been 
established. 

4 .  Regarding  LIDAR  Shipboard  Implementation 

•  Wavelengths  longer  than  1.4  ym  provide  a  reduc¬ 
tion  of  over  500  in  eye  safety  hazards  for  the 
same  laser  pulse  energy. 

•  Increasing  the  transmitted  beam  diameter  to 
30.5  cm  (12  inches)  from  0.8  cm  (ASL  Model-1 
visioceilometer)  reduces  the  transmitted  beam 
energy  density  at  the  transmitter  by  a  factor 
of  1450. 

•  Army  AN/TVQ-2  (GLLD)  and  Navy  A-6  TRAM  laser 
designator  systems  are  promising  candidates  as 
sources  for  reliable,  high  performance,  MIL- 
SPEC  hardware  modules. 

•  A  prototype  system  using  a  wavelength  shifted 
1.06  ym  laser  (shifted  to  1.54  ym  using  a  sti¬ 
mulated  Raman  scattering  cell  in  tandem  with 
the  1.06  ym  source)  and  a  transmit/receive  op¬ 
tics  diameter  of  30.5  cm  should  provide  good 
performance  for  ranges  in  excess  of  5  km  while 
remaining  eye-safe. 

•  A  two-color  LIDAR  system  should  be  developed. 
Simultaneous  near  IR  (1.54  ym)  and  10.6  ym 
shipboard  measurements  will  be  required  to: 

-  generate  needed  information  to  supplement 
the  inadequate  exsisting  marine  aerosol 
data  base 

-  evaluate  and  develop  a  system  required  for 
reliable  routine  shipboard  measurement  of 
aerosol  extinction  in  the  3-5  ym  and  8-12  ym 
FLIR  operating  bands. 
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2 

CHARACTERIZATION  OF  THE  MARINE  ENVIRONMENT 


2.1  MARINE  AEROSOL  MODELS 

In  this  section  we  will  be  interested  in  the  relation¬ 
ship  of  aerosol  attenuation  at  infrared  wavelengths  to  that 
in  the  visible  or  photopic  region  which  can  be  derived  from 
a  simple  visibility  measurement. 

Aerosol  attenuation  at  optical  wavelengths  (visible 
light)  can  be  directly  related  to  the  meteorological  visi¬ 
bility  or  visual  range  through  Koschmeider ' s  law: 


where  a .  is  the  aerosol  extinction  coefficient  and 
A 

the  visual  range  or  visibility.  This  range  corresponds  to 
the  distance  over  which  an  observer  can  barely  discern  a 
distant  object.  The  contrast  between  the  object  and  the 
surrounding  background  in  the  observer's  f ield-of-view  is 

taken  to  be  0.02  by  convention.  When  R^  is  expressed  in 

.  -1  v 
km  then  aA  has  units  of  km 

The  simplest  relationship  or  model  between  visible 
and  infrared  aerosol  extinction  which  has  been  used  in  the 
past  has  the  form: 


aiR  =  °VIS  X 


-n 


where  oIR  and  ®VIS  are  aerosol  extinction  coefficients  for 
infrared  and  visible  wavelengths,  respectively,  X  is  wave¬ 
length  and  n  is  an  exponent  slightly  greater  than  unity. 
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This  type  of  scaling  between  visible  and  infrared 
aerosol  extinction  has  been  shown  to  be  better  suited  for 
describing  extinction  by  continental  rather  than  marine 
aerosols  where  a  larger  variation  in  the  ratio  of  visible 
to  infrared  extinction  is  seen  than  for  the  continental 
case.  In  an  effort  to  better  quantify  the  variations  both 
in  magnitude  of  visible  aerosol  extinction,  i.e.,  visi¬ 
bility,  and  also  the  ratio  of  infrared  to  visible  extinc¬ 
tion,  models  have  been  developed  for  marine  aerosol  scat¬ 
tering  which  relate  the  large  variations  in  observed  ex¬ 
tinction  directly  to  changes  in  the  shape  of  the  marine 
aerosol  distribution.  The  modeled  marine  aerosol  distribu¬ 
tion  is  then  used  in  a  Mie  scattering  calculation  to  predict 
the  magnitude  of  aerosol  extinction  at  the  wavelength  of 
interest.  Not  only  does  the  density  of  aerosol  particles 
per  unit  volume  scale  with  visibility  (visible  extinction) , 
which  is  predominately  determined  by  particle  sizes  with 
radii  of  a  few  tenths  of  a  micron,  but  the  shape  of  the 
distribution  also  changes,  depending  upon  additional  input 
parameters  such  as  windspeed  and  relative  humidity.  Since 
particles  with  radii  comparable  to  the  wavelength  under 
consideration  are  most  effective  in  scattering  that  wave¬ 
length,  the  decreased  fall-off  of  particle  number  density 
with  increasing  particle  size  in  the  marine  aerosol  model 
results  in  infrared  scattering  more  nearly  comparable  with 
that  in  the  visible.  In  this  regard  recent  marine  aerosol 
models  are  more  consistent  with  both  maritime  optical  and 
infrared  transmission  measurements  as  well  as  aerosol  dis¬ 
tribution  measurements  performed  in  maritime  locations. 

The  development  of  marine  aerosol  models  has  been  the 
subject  of  much  interest  in  recent  years  and  has  been  dis¬ 
cussed  recently  in  detail  by  Richter  and  Hughes  [4]  and  by 
Katz  15] .  Only  the  essential  features  of  the  models  will 
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be  described  here  since  we  are  primarily  interested  in  the 
predicted  ranges  of  infrared  and  visible  aerosol  scattering 
values  obtained  using  these  models  and  not  in  the  develop¬ 
ment  or  detailed  validation  of  the  models  themselves.  We 
will  also  summarize  some  recent  infrared  and  visible  aerosol 
scattering  measurements  so  that  together  with  the  range  of 
predicted  scattering  values  obtained  with  the  models,  we 
are  likely  to  bracket  a  range  of  values  representative  of 
the  marine  environment. 

The  comparisons  of  measured  and  predicted  aerosol 
scattering  for  the  marine  environment  ti.  be  discussed  later 
should  by  no  means  be  taken  to  be  exhaustive  or  conclusive. 
Within  the  scope  of  the  effort  presented  here  such  compari¬ 
sons  are  rather  intended  to  contrast  certain  features  of 
recent  modeling  developments  and  to  review  a  few  selected 
examples  of  marine  aerosol  scattering  measurements  to  pro¬ 
vide  a  reasonable  assessment  of  the  environment  in  which  a 
shipboard  LIDAR  system  will  operate. 

Aerosol  models  for  the  marine  environment  are  charac¬ 
terized  by  an  admixture  of  two  components  (continental  and 
maritime)  .  The  work  of  Bamhardt  and  Streete  [6]  ,  the 
earlier  version  of  LOWTRAN  3B  [7]  and  the  recent  LOWTRAN  5 [8] 
model  are  all  based  on  this  combination  of  two  components/ 
but  various  descriptions  are  used  for  each  in  the  different 
models.  See  reference  4  for  a  discussion  of  these  differ¬ 
ences. 

The  LOWTRAN  models  use  a  mixture  of  oceanic  and  con¬ 
tinental  components  in  a  fixed  ratio  with  the  total  number 
density  normalized  by  visibility.  The  relative  contribu¬ 
tions  of  each  component  and  their  mathematical  form  differ 
in  the  two  versions. 


11 


The  model  of  Wells,  Gal  and  Munn  (WGM)  [9]  and  the 
Katz/Ruhnke  (K/R)  model  [10]  both  contain  continental  and 
oceanic  components  as  in  the  LOWTRAN  models  but  not  in  fix¬ 
ed  proportions  to  each  other.  In  the  WGM  and  K/R  models 
the  relative  contribution  of  the  oceanic  component  depends 
on  windspeed  and  relative  humidity.  The  major  distinction 
between  the  WGM  and  K/R  models  has  to  do  with  a  difference 
in  altitude  dependency,  which  is  an  exponential  particle 
density  decrease  in  the  WGM  model  and  a  particle  size  de¬ 
pendent  form  in  the  K/R  model  that  provides  a  more  uniform 
height  distribution  of  larger  particles  during  high  wind 
conditions. 

In  the  K/R  model  the  continental  to  maritime  mixing 
ratio  is  determined  by  the  saturation  ratio  (relative 
humidity  divided  by  100)  whereas  it  is  an  adjustable  input 
parameter  in  the  WGM  model.  The  particle  growth  factor 
with  relative  humidity  was  changed  in  the  K/R  model  from 
that  of  Barnhardt  and  Streete  (also  used  in  the  WGM  model) 
to  that  devised  by  Fitzgerald  [11] .  In  the  WGM  model  the 
calculated  aerosol  size  distribution  is  normalized  by  the 
observed  visual  range  but  this  feature  was  removed  in  the 
K/R  model  due  to  the  questionable  precision  and  utility 
of  shipboard  meteorological  observer  visibility  measurements. 
Results  obtained  with  the  K/R  model  will  be  compared  to 
those  obtained  with  the  two  versions  of  LOWTRAN  in  the  sub¬ 
sections  that  follow. 

The  recently  developed  K/R  model  is  an  outgrowth  of 
the  WGM  approach  and  in  addition  to  the  differences  cited 
above  it  incorporates  changes  in  the  values  of  some  of  the 
constants  used  in  the  WGM  particle  size  distribution  for¬ 
mulae.  Consequently  the  more  recent  K/R  model  rather  than 
the  earlier  WGM  approach  will  be  considered  in  comparisons 
to  the  LOWTRAN  models  and  to  experimental  data. 
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A  comparison  of  the  aerosol  particulate  distributions 
calculated  separately  with  the  oceanic  and  continental  com¬ 
ponents  of  the  LOWTRAN  3B  model  and  their  sum  (the  LOWTRAN 
3B  maritime  distribution)  is  shown  in  Figure  1.  As  one  can 
see  in  the  figure  the  continental  component  dominates  the 
LOWTRAN  maritime  distribution  for  particle  sizes  below 
about  0.5  urn  while  the  oceanic  component  is  larger  for 
sizes  greater  than  this  value.  The  K/R  model  on  the  other 
hand,  predicts  a  decreased  concentration  of  laraer  par¬ 
ticles  and  consequently  a  smaller  ratio  of  aTO/o,.Te  (ratio 
of  infrared  to  visible  aerosol  extinction) ,  than  does  LOW¬ 
TRAN  for  3.5  m/s  windspeed,  growing  to  an  increased  value 
compared  to  LOWTRAN  for  a  20  m/s  windspeed.  The  two  model 
distributions  are  seen  to  coincide  for  an  intermediate 
windspeed  value  around  10  m/s. 


2.1.1  WINDSPEED  DEPENDENCE 


The  impact  of  the  similarities  and  distinctions  be¬ 
tween  the  LOWTRAN  models  and  the  K/R  approach  can  be  readily 
seen  by  examining  the  ratio  of  infrared  aerosol  extinction 
to  visible  extinction  predicted  by  each  model  for  different 
values  of  relative  humidity  and  windspeed.  Figure  2  shows 
data  taken  from  reference  5  plotted  as  the  ratio  of  infra¬ 
red  to  visible  aerosol  scattering  coefficient,  aIR/aVIS 
versus  wind  speed  in  m/s. 


Figure  2(A)  shows  a  comparison  of  the  ratio  a3_5/aVIg 
for  the  3-5  ym  IR  band  for  the  three  models  previously  dis¬ 
cussed,  namely  LOWTRAN  3B,  LOWTRAN  5  and  the  Katz/Rhunke 
(K/R)  model.  The  LOWTRAN  3B  calculation  is  independent  of 
variations  in  windspeed  and  relative  humidity  and  predicts 
a  constant  value  for  aIR/aVIg  for  the  3-5  ym  band  (Figure 
2(A))  and  a  smaller  constant  value  for  the  8-12  ym  band 
(Figure  2(B)). 
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PARTICLE  SIZE  C  /j  m) 


FIGURE  1.  AEROSOL  DISTRIBUTION  FOR  LOWTRAN  AND  KATZ- 
RUHNKE  MODELS . 
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FIGURE  2.  RATIO  OF  INFRARED  TO  VISIBLE  AEROSOL 
EXTINCTION  VERSUS  WINDSPEED  FOR  LOW- 
TRAN  AND  K/R  MARITIME  AEROSOL  MODELS: 
(A)  3-5  ym  REGION;  (B)  8-12  um  REGION. 
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The  LOWTRAN  5  model  predicts  a  significant  variation 
of  25%  in  aIR/aVIS  for  a  chani?e  of  RH  from  80%  to  95%  for 
the  3-5  uni  band  and  about  the  same  relative  change  for  the 
8-12  ym  band  although  the  actual  values  of  aT_/a,.TC  for 

IK  V  Xb 

8-12  ym  are  only  about  a  third  of  those  for  the  3-5  ym  band. 

The  K/R  model  by  contrast  exhibits  a  smaller  variation 
in  ajR/aVIS  RH  than  LOWTRAN  5  (only  5%  to  12%  over  the 

range  of  wind  speeds  from  5  m/s  to  15  m/s  for  the  3-5  ym 
band) .  The  absolute  values  of  aIR/ayIS  predicted  by  the 
K/R  model  only  approach  the  LOWTRAN  5  predictions  for  wind- 
speeds  in  excess  of  11  to  12  m/s  and  always  fall  below  the 
95%  RH  LOWTRAN  5  calculation  for  the  3-5  ym  band  while 
matching  the  8-12  ym  value  for  13  m/s  windspeed  and  exceed¬ 
ing  it  for  higher  values.  The  region  of  overlap  of  pre¬ 
dictions  for  the  two  models  is  seen  to  be  confined  to  higher 
windspeeds  >  11  m/s  and  the  K/R  model  predictions  for  a  / 

X  K 

aVIS  are  always  lower  for  windspeeds  less  than  this  value. 

2.1.2  RELATIVE  HUMIDITY  DEPENDENCE 

The  sensitivity  of  infrared  aerosol  scattering  to 
changes  in  relative  humidity  for  the  three  marine  aerosol 
models  can  be  examined  by  again  plotting  the  ratio  o  D/a..TO, 
but  in  this  case  versus  relative  humidity  (RH) .  The  values 
of  aIR/CTVIS  predicted  by  each  model  overlap  for  the  range 
of  RH  values  between  70%  and  90%  and  for  the  range  of  wind- 
speeds  between  10  m/s  and  15  m/s.  As  shown  in  Figure  3, 
the  predictions  obtained  with  LOWTRAN  3B  are  independent 
of  changes  in  RH.  LOWTRAN  5  predictions  on  the  other  hand 
bracket  the  LOWTRAN  3B  values,  predicting  lower  values  below 
about  78%  RH  for  the  3-5  ym  case  and  below  85%  RH  for  8-12 
ym,  and  higher  values  of  ajR/avis  each  case  for  higher 
values  of  RH.  The  K/R  model  predictions  by  comparison  span 
the  range  of  LOWTRAN  5  predictions  for  each  wavelength 
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FIGURE  3.  RATIO  OF  INFRARED  TO  VISIBLE 

AEROSOL  EXTINCTION  VERSUS  RELA¬ 
TIVE  HUMIDITY  FOR  LOWTRAN  AND 
K/R  MARITIME  AEROSOL  MODELS: 

(A)  3-5  yin  REGION;  (B)  8-12  ym 
REGION. 
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region  (3-5  urn  shown  in  Figure  3(A)  and  8-12  urn  shown  in 
Figure  3(B))  for  windspeed  values  between  10  m/s  and  15  m/s. 
In  both  wavelength  regions  the  K/R  model  predictions  for 
5  m/s  windspeeds  are  substantially  lower  than  either  LOWTRAN 
model  prediction. 

Within  the  combined  range  of  10  m/s  to  15  m/s  wind- 
speed  and  70%  to  95%  RH  the  models  are  thus  seen  to  pre¬ 
dict  comparable  results.  Values  of  aTn/a..__  predicted  bv 
LOWTRAN  3B  are  constant,  independent  of  both  windspeed  and 
RH.  LOWTRAN  5  is  one  step  more  sophisticated  by  incorporat¬ 
ing  the  RH  dependence  seen  in  Figure  3,  but  does  not  include 
a  component  representing  a  windspeed  dependent  sea-spray 
contribution  as  does  the  K/R  model.  As  is  evident  in  Fi¬ 
gures  2  and  3,  the  K/R  model  incorporates  both  a  windspeed 
and  RH  dependence  and  therefore  should  be  more  representa¬ 
tive  of  observations  by  allowing  for  independent  variations 
in  both  of  these  parameters.  The  values  of  aT.T„/aT_,  pre- 
dieted  by  the  K/R  model  as  always  lower  than  the  LOWTRAN 
models  for  windspeeds  less  than  11  to  12  m/s  due  to  the 
predominance  of  the  continental  aerosol  component  in  the 
model  for  low  windspeeds. 


2.1.3  WAVELENGTH  DEPENDENCE 


Based  on  the  discussions  presented  in  the  previous  two 
sections  the  region  of  comparability  of  results  obtained 
with  the  three  aerosol  models  under  consideration  is  con¬ 
fined  to  the  combined  ranges  of  10-15  m/s  windspeeds  and 
70%-90%  RH.  Within  this  region,  the  values  of  aij*/aVis  ^or 
the  3-5  urn  predicted  by  the  three  models  vary  by  a  little 
more  than  a  factor  of  2  from  0.38  to  0.78.  The  same  type 
of  comparison  shows  about  a  factor  of  3  variation  for  the 
8-12  um  region,  but  lower  overall  values  of  the  ratio  ag_^2// 

The  values  of  oIR/°VIS  predicted 


aVIS  between  0.12  and  0.38. 
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by  LOWTRAN  3B  are  fixed  at  0.55  and  0.21  for  3-5  ym  and 
8-12  ym  respectively,  independent  of  both  windspeed  and  RH. 
LOWTRAN  5  shows  the  variation  with  RH  seen  in  Figure  3  but 
is  independent  of  windspeed  while  the  K/R  model  contains 
both  dependencies  as  previously  noted. 


The  wavelength  dependence  of  the  K/R  model  and  its ' 
dependence  upon  windspeed  is  demonstrated  in  Figure  4 .  For 
low  windspeeds  of  5  m/s  the  extinction  ratio  aT_/a.7_c,  is 
40%  at  4  ym  of  its'  value  at  1.06  ym  and  only  15%  of  this 
value  at  10  ym.  For  increased  windspeeds  the  magnitude  of 
aIR/aVIS  increases  and  the  decrease  in  this  ratio  for  longer 
wavelengths  becomes  weaker.  This  behavior  reflects  the 
relative  growth  of  the  wind-driven  sea  spray  component  in 
the  K/R  aerosol  model,  causing  aerosol  extinction  in  the 
infrared  to  increase  relative  to  visible  scattering  as 
greater  numbers  of  larger  sea-spray-derived  particles  are 
generated.  The  graph  of  aIR/aVIS  tends  to  flatten  out  with 
increasing  wavelengths  just  as  the  actual  particle  distri¬ 
butions  shown  in  Figure  1  show  a  slower  fall-off  of  dN/dR 
for  larger  particle  sizes. 

The  aerosol  model  predictions  discussed  in  this  section 
will  be  compared  to  experimental  data  in  Section  2.2. 


2.1.4  CLIMATOLOGIES 

For  our  purposes  it  is  useful  to  consider  results  ob¬ 
tained  from  statistical  analyses  of  maritime  meteorological 
data  collected  during  periods  of  several  years  which  have 
been  used  in  developing  E-0  systems  performance  estimates. 
Exsisting  weather  ship  data  have  been  used  by  Davis  et  al 
[12],  Katz  et  al  [13]  and  recently  by  Katz  and  Goroch  [14] 
to  produce  climatological  statistics.  These  have  been  com¬ 
bined  with  atmospheric  effects  models  to  generate  statisti¬ 
cal  performance  estimates  for  E-0  systems.  Typically 
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studies  such  as  these  have  concentrated  on  performance 
estimates  for  Forward  Looking  Infrared  (FLIR)  systems 
operating  in  the  3-5  um  and  8-12  ym  regions.  Aerosol  ex¬ 
tinction  models  such  as  those  in  LOWTRAN  and  the  K/R  model 
as  well  as  molecular  absorption  calculations  are  used  in 
relating  system  performance  to  the  weather  data. 

To  the  extent  that  "climatological-like"  statistics  of 
aerosol  extinction  at  infrared  wavelengths  are  available  as 
published  intermediate  results  in  studies  such  as  references 
12-14,  they  can  be  used  to  assess  the  range  of  aerosol 
scattering  values  to  be  expected  at  different  wavelengths 
and  for  different  maritime  locations. 

Such  "climatological"  results  are  only  as  valid  as  the 
specific  aerosol  extinction  model  used.  There  is  evidence 
that  factors  related  to  air  mass  character  and  length  of 
long-term  averaging  time  for  windspeed,  which  are  not 
explicitly  contained  in  the  aerosol  models  are  correlated 
with  observed  extinction  values  and/or  measured  maritime 
aerosol  distributions. 

2.2  MARINE  AEROSOL  MEASUREMENTS 

2.2.1  LONG-PATH  OPTICAL  TRANSMISSION  MEASUREMENTS 

Several  measurement  programs  have  been  undertaken 
at  various  maritime  locations  throughout  the  world  during 
recent  years.  Extended  over-water  measurement  programs 
were  begun  around  1975  by  Physical  Laboratory  TNO  in 
Holland  [15] ,  more  recently  by  the  Electronics  Research 
Laboratory,  Defense  Science  and  Technology  Organization, 
South  Australia  [16-18]  and  by  the  Pacific  Missile  Test 
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Center  (PMTC)  at  San  Nicolas  Island  (SNI),  California  [19]. 
Representative  results  obtained  from  these  measurement  pro¬ 
grams  were  recently  compared  in  the  context  of  a  Workshop  on 
Atmospheric  Water  Vapor  [20].  In  many  of  the  cases  examined 
in  reference  20  comparisons  of  the  optical  measurements  to 
models  such  as  LOWTRAN  pointed  out  probable  deficiencies 
in  the  atmospheric  molecular  absorption  (water  vapor) 
modeling  rather  than  the  aerosol  scattering  components  of 
the  model.  Most  of  the  data  considered  were  measured  under 
relatively  high  visibility  conditions,  consequently  the 
aerosol  attenuations  were  correspondingly  low. 

Not  considered  in  reference  20  but  relevant  to  the 
present  considerations  are  selected  data  collected  under 
the  OPAQUE  program  [21].  In  particular,  data  collected  at 
the  Netherlands  OPAQUE  site  near  the  North  Sea  are  often 
strongly  subjected  to  the  influence  of  maritime  conditions 
arising  during  on-shore  wind  conditions.  These  data  have 
been  used  in  the  development  of  the  Electro-Optical  Systems 
Atmospheric  Effects  Library  (EOSAEL)  [22].  Selected  ex¬ 
amples  from  this  analysis  will  be  discussed  below. 

Long-Path  infrared  transmission  measurements  have  also 
been  performed  by  the  Naval  Research  Laboratory  at  a  variety 
of  locations  which  have  included  coastal  sites  at  Cape  Ca¬ 
naveral  Air  Force  Station,  Florida  and  San  Nicolas  Island, 
California  [23-24] .  Representative  results  from  these 
experiments  taken  during  on-shore  wind  conditions  providing 
maritime  air  along  the  transmission  measurement  path  will 
also  be  reviewed. 

Recent  Marine  Aerosol  Distribution  Measurements  will 
be  discussed  in  Section  2.2.2. 
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2. 2. 1.1  Physical  Laboratory  TNO  (Holland)  Long-Path 
Infrared  Transmission  Measurements 

The  Physical  Laboratory  TNO,  Netherlands,  has  been 
conducting  maritime  atmospheric  transmission  measurements 
at  two  locations  on  the  North  Sea  coast  [15] .  One  is 
across  an  estuary  of  the  Scheldt  near  Flushing,  along  a 
6  km  path.  More  recently,  an  18.2  km  path  between  a  shore 
site  and  an  artificial  island  in  the  North  Sea  has  been 
used.  Two  bands  were  studied  in  the  first  set  of  measure¬ 
ments  conducted  during  1974  and  1975,  namely,  3.4  ym  to 
4.1  ym  and  8.2  ym  to  13.8  ym.  Uncooled  pyroelectric  de¬ 
tectors  were  used.  More  recently  in  the  longer  path 
measurements ,  cooled  detectors  were  employed  and  several 
individual  filter  bands  in  each  of  the  two  regions  were 
studied. 

Figures  5  and  6  show  transmission  data  for  these  ex¬ 
periments  plotted  against  atmospheric  water  vapor  density 
along  with  LOWTRAN  calculations  for  the  18.2  km  path.  The 
data  for  the  mid-IR  region  between  3.4  ym  and  4.1  ym 
shown  in  Figure  5  display  a  substantially  smaller  decrease 
in  transmission  with  increasing  water  vapor  density  than  do 
corresponding  data  for  8  ym  and  13  ym  shown  in  Figure  6. 

It  is  interesting  to  note  that  the  ranges  of  transmission 
values  are  not  significantly  different  for  the  two  bands 
although  the  slopes  of  the  lines  fitted  to  the  data  are. 

The  measured  data  for  the  long  wavelength  band  shown 
in  Figure  6  show  greater  transmission  than  the  LOWTRAN 
model  predicts  for  the  18.2  km  path.  The  data  for  the 
mid-IR  band  shown  in  Figure  5  show  greater  scatter  but 
appear  to  be  in  better  agreement  with  the  LOWTRAN  model 
for  the  band  between  3.4  ym  and  4.1  ym. 
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FIGURE  5.  ATMOSPHERIC  TRANSMISSION  IN  THE  3.4  pm  TO  4.1  pm  BAND  MEASURED 
OVER  18.2  km  AND  6  km  PATHS  (TNO-NETHERLANDS ) . 


FIGURE  6.  ATMOSPHERIC  TRANSMISSION  IN  THE  8  um  TO  13  pm  BAND  MEASURED  OVER 
18.2  km  AND  6  km  PATHS  (TNO-NETHERLANDS) . 


These  data  were  taken  for  conditions  of  moderate 
visibility,  i.e.  R^  >  15  km  (aVIS  <  0.26  km  1).  The  rea¬ 
sonably  good  agreement  of  the  data  with  the  LOWTRAN  3B  cal¬ 
culations  for  Ry  =  16  km  shown  for  comparison  in  Figure  5 
indicates  that  the  sum  of  contributions  due  to  molecular 
absorption  and  aerosol  extinction  are  being  accurately  pre¬ 
dicted  by  the  model.  The  variability  of  the  data  for  the 
18.2  km  path  measurements  for  the  range  of  atmospheric 
water  vapor  density  around  75-80  precipitable  mm  per  18.2 
km  may  well  be  accounted  for  by  variability  in  the  actual 
aerosol  scattering  not  modeled  by  LOWTRAN  3B  but  due  to 
wind  and  relative  humidity  variations  which  are  taken  into 
account  in  the  WGM  and  K/R  models.  For  the  low  values  of 
absolute  humidity  corresponding  to  these  measurements  (i.e. 

3 

3  to  6  gm/m  of  water  vapor)  the  coefficient  for  molecular 
absorption  should  be  only  about  10%  of  that  for  aerosol 
extinction.  Accordingly  the  general  agreement  between  the 
data  shown  in  Figure  5  and  the  LOWTRAN  3B  prediction  is 
primarily  determined  by  agreement  between  the  measurements 
and  the  LOWTRAN  3B  aerosol  model  predictions. 

The  TNO  data  shown  in  Figure  6  for  the  8-13  ym  are 
not  as  useful  as  the  mid-infrared-band  data  in  verifying 
the  LOWTRAN  aerosol  model  predictions.  Here  the  relative 
contributions  of  molecular  absorption  and  aerosol  scatter¬ 
ing  are  more  nearly  comparable,  the  former  being  about  a 
factor  of  six  larger  than  for  the  data  shown  in  Figure  5 
while  aerosol  scattering  in  the  8-13  ym  band  should  be  less 
than  30%  of  that  in  the  shorter  wavelength  band  for  the 
atmospheric  conditions  corresponding  to  most  of  the  data 
shown  in  Figures  5  and  6,  (i.e.  15-18  km  visibilities,  60% 
to  80%  relative  humidities  and  wind  speeds  >  10  m/s)  [25]. 

For  the  conditions  corresponding  to  most  of  the  TNO 
measurements,  the  LOWTRAN  3B,  LOWTRAN  5  and  K/R  aerosol 
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models  will  give  roughly  comparable  predictions.  Differences 
in  predictions  by  the  various  models  only  become  apparent 
for  windspeeds  below  12  m/s  (see  Figure  2)  where  the  K/R 
model  predicts  smaller  values  of  aTD/o,7TO  than  does  either 
version  of  LOWTRAN,  and/or  for  relative  humidities  <  70% 
or  >  90%. 


2. 2. 1.2  Australian  Defense  Science  Establishment 

Long-Path  Infrared  Transmission  Measurements 

The  Electronics  Research  Laboratory  of  the  Australian 
Defense  Scientific  Technical  Office  carried  out  a  measure¬ 
ment  program  using  maritime  paths  located  at  Victor  Harbor, 
80  km  south  of  Adelaide  in  1977  and  1978.  The  path  lengths 
used  were  5.03  km  and  9.05  km.  A  large  reentrant  black- 
body  cavity  at  a  temperature  of  980  K  was  used  as  a  source 
and  combined  with  an  all-reflective  receiver  by  using 
three  detectors  (PbS,  InSb,  and  HgCdTe)  and  various  filters 
to  cover  selected  spectral  bands  between  1  pm  and  13  um 
U8]  . 

Figures  7  and  8  show  samples  of  data  collected  for 
each  of  the  two  paths  during  this  program.  Here  transmis¬ 
sion  in  the  band  between  4.41  pm  and  5.4  ym  is  plotted 
against  water  vapor  concentration  in  gm/m3 .  The  data  were 
collected  during  very  high  visibility  conditions  for  which 
the  visual  range  exceeded  96  km.  The  figures  show  compari¬ 
sons  of  LOWTRAN  model  calculations  with  the  measured  trans¬ 
mission  values.  The  experimental  data  are  at  variance  with 
the  LOWTRAN  calculation  for  water  vapor  densities  greater 
than  about  9  gm/m3.  Figures  9  and  10  show  comparable  data 
for  the  8.2  ym  to  11.8  ym  band  for  the  same  two  paths, 
again  compared  with  the  LOWTRAN  model.  Here  the  agreement 
between  measurement  and  calculation  for  the  shorter  path  is 
better  than  that  for  the  shorter  wavelength  band  shown  in 
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FIGURE  9.  ATMOSPHERIC  TRANSMISSION  IN  THE  8.2  vim  TO  11.8  pm  BAND 
MEASURED  OVER  A  9.05  km  PATH  (DSTO-AUSTRALIA) . 
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FIGURE  10.  ATMOSPHERIC  TRANSMISSION  IN  THE  8.2  pm  TO  11.8  pm  BAND 
MEASURED  OVER  A  5.03  km  PATH  (DSTO-AUSTRALIA). 


Figures  7  and  8.  However,  for  the  9.05  km  path  LOWTRAN 
significantly  underestimates  the  transmission. 

The  data  shown  in  Figures  7  through  10  are  of  interest 
in  testing  the  validity  of  the  molecular  absorption  compo¬ 
nents  of  LOWTRAN  but  are  not  particularly  useful  in  verify¬ 
ing  aerosol  models  since  the  actual  aerosol  extinction  com¬ 
ponent  of  the  measurements  is  vanishingly  small  (R^  >  96  km 
in  all  cases) .  All  outdoor  transmission  measurements  will 
contain  contributions  from  each  source  of  attenuation. 
Hopefully,  selection  of  certain  sub-sets  of  data  can  be 
made  depending  upon  which  components  of  a  transmission 
model  are  of  interest.  The  Australian  DSTO  data  are  useful 
for  molecular  absorption  model  validation  but  not  for  the 
topic  of  present  concern  in  this  work,  i.e.,  aerosol  model 
comparisons.  The  TNO  data  on  the  other  hand  which  are 
shown  in  Figures  5  and  6  were  collected  under  moderate 
visibility  <RV  *\»  15  km)  and  low  water  vapor  (3-6  torr  H20) 
conditions  and  are  therefore  more  useful  in  aerosol  model 
comparisons,  particularly  the  data  for  the  3-5  ym  region. 
Data  in  which  the  relative  contribution  of  molecular  ab¬ 
sorption  and  aerosol  extinction  are  comparable,  are  typi¬ 
cally  less  useful  for  model  comparison  purposes,  except  in 
cases  where  the  molecular  absorption  coefficients  are  very 
well  known.  The  residual  attenuation  is  then  a  reliable 
measure  of  the  aerosol  contribution  to  total  attenuation. 

2.2.1. 3  OPAQUE  Measurement  Program  Long-Path  Infrared 
Transmission  Measurements 

An  extensive  measurement  program  at  several  sites  in 
Western  Europe  was  begun  under  the  OPAQUE  program  in  1976 
and  continued  into  1980  [21] .  Selected  data  from  the 
Netherlands  measurement  site  located  at  the  Hague  near  the 
North  Sea  coast  were  analyzed  in  reference  22.  The  data 
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available  for  that  study  were  collected  during  the  period 
from  1  March  1977  to  31  May  1977.  At  this  time  of  year  the 
Netherlands  is  subject  to  a  wide  range  of  weather  conditions 
varying  from  cold  temperatures  and  dry  conditions  resulting 
from  polar  airmasses  to  mild,  moist  conditions  associated 
with  maritime  air  masses. 

Examples  of  the  OPAQUE  data  analyzed  in  reference  22 
are  shown  in  Figures  11  and  12  where  they  are  compared 
to  LOWTRAN  predictions  and  to  a  water  haze  model  whose  de¬ 
velopment  is  described  in  reference  22.  Supporting  meteo¬ 
rological  data  for  the  optical  measurements  shown  here  are 
not  readily  available  without  recourse  to  magnetic  data 
tapes  obtainable  from  the  OPAQUE  program.  Detailed  analy¬ 
sis  of  these  data  and  comparison  to  LOWTRAN  5  and  K/R  model 
predictions  is  beyond  the  scope  of  the  present  effort, 
however  the  trends  apparent  in  the  data  shown  in  Figures 
11  and  12  and  the  comparisons  to  the  maritime  and  rural 
components  of  LOWTRAN  3B  are  instructive  to  note.  The  data 
for  the  ratio  a3_5/avis  ^etween  the  LOWTRAN  calculations 
for  the  two  types  of  atmospheres  and  are  in  fair  agreement 
with  the  water  haze  model  described  in  Reference  22.  A  large 
amount  of  scatter  is  evident  in  the  data  however.  The  data 
for  the  8-12  urn  region  shown  in  Figure  12  show  quite  a 
different  trend  indicating  that  as-i2/,aviS  is  9enerally 
larger  than  the  predictions  obtained  with  either  version 
of  LOWTRAN  or  the  SAI  water  haze  model.  The  ranges  of 
values  for  anc*  a8-12  s^own  in  't^ie  figures  are  still 

useful  in  assessing  the  range  of  scattering  values  at  these 
wavelengths  to  be  expected  for  maritime  conditions. 
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flER  EXT  COEF  6  3-5 


FIGURE  11.  COMPARISON  OF  VARIOUS  AEROSOL  MODELS  WITH 
OPAQUE  DATA  AEROSOL  EXTINCTION  COEFFICIENT 
§  0.55  um  VERSUS  AEROSOL  EXTINCTION  COEFFI¬ 
CIENT  @  3-5  um  DERIVED  FROM  OPAQUE  NETHER¬ 
LANDS  3/77-5/77  DATA. 
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PER  EXT  COEF  8  8-12 


FIGURE  12.  COMPARISON  OF  VARIOUS  AEROSOL  MODELS  WITH 
OPAQUE  DATA  AEROSOL  EXTINCTION  COEFFICIENT 
@  0.55  yrn  VERSUS  AEROSOL  EXTINCTION  COEFFI¬ 
CIENT  @  8-12  ym  DERIVED  FROM  OPAQUE  NETHER¬ 
LANDS  3/77-5/77  DATA. 
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2. 2. 1.4  Pacific  Missile  Test  Center  (PMTC) /Naval 
Weapons  Center  Optical  Signatures  Program 
(NWC-OSP)  Long-Path  Infrared  Transmission 
Measurements 

Overwater  long-path  infrared  transmission  measurements 

using  2.4  km  and  4.07  km  paths  have  been  performed  at  San 

Nicolas  Island,  California  starting  in  1978  [26] .  A  small 

quantity  of  preliminary  data  generated  by  this  program  was 

analyzed  and  described  in  reference  22.  Examples  of  volume 

extinction  coefficients  at  1.06  pm,  3.7  -  3.9  pm  and  4.8  - 

5.0  pm  plotted  versus  0.55  pm  extinction  coefficients  from 

this  limited  set  of  data  are  shown  in  Figures  13,  14  and 

15  respectively.  The  reasonably  good  linear  curve  fits  to 

the  aTT3/o.  c  ratios  shown  in  these  figures  indicate  that  a 
IK  VXS  ^ 

simple  approximate  relationship  of  the  form  aIR  v  aViS 
might  be  used  to  relate  infrared  to  visible  aerosol  attenua¬ 
tion.  When  much  larger  samples  of  data  are  considered 
greater  scatter  is  typically  observed  and  a  power  law 
relationship  no  longer  fits  the  data  very  well.  It  is 
likely  that  this  observation  corresponds  to  the  super- po¬ 
sition  of  data  sets  corresponding  to  different  values  of 
n  in  a  power  law  relationship  which  may  be  related  to 
different  air  masses.  That  is,  a  given  value  of  n  may  pro¬ 
vide  a  good  fit  to  data  such  as  shown  in  Figures  13-15  under 
the  influence  of  a  constant  air  mass,  whereas  a  different 
value  of  n  may  be  appropriate  for  different  air  masses. 

When  the  data  for  several  air  masses  are  superimposed  on  a 
single  plot  of  crIR  versus  aVIg,  no  single  value  of  n  fits 
the  data  and  the  more  sophisticated  aerosol  models  (e.g. 

K/R  model  dependence  on  windspeed  and  RH)  are  required. 
Recently  a  more  comprehensive  statistical  analysis  of 
data  generated  under  the  PMTC/OSP  measurement  program  be¬ 
tween  1978  and  1979  has  been  performed  by  Katz  and  DeBold 
[27] .  In  this  analysis  data  were  compared  to  the  K/R, 
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TAKEN  AT  SAN  NICOLAS  ISLAND  BETWEEN  4/27/78  AND  5/18/78. 
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FIGURE  15.  VOLUME  EXTINCTION  COEFFICIENTS  AT  4. 5-5.0  yin  AND  0.55  ym  FOR  OSP  DATA  TAKEN 
AT  SAN  NICOLAS  ISLAND  BETWEEN  4/27/78  AND  5/18/78. 


LOWTRAN  3B  Maritime  and  GAP  (low  visibility)  [22]  models. 
Figure  16  shows  an  example  of  comparisons  of  the  OSP-SNI 
data  for  the  ratios  aIR/ayjg  for  different  wavelength  bands 
plotted  versus  windspeed.  The  extinction  data  were  divided 
into  windspeed  intervals  of  2  m/s.  The  error  bars  indi¬ 
cate  the  standard  deviations  about  the  mean  value  for  each 
windspeed  group.  The  LOWTRAN  3B  model  predictions  are  in 
agreement  with  the  experimental  data  (solid  lines)  for 
windspeeds  near  10  m/s.  The  broken  lines  show  calculations 
performed  with  the  K/R  model  which  are  seen  to  be  in 
reasonable  agreement  with  the  experimental  data.  The 
LOWTRAN  calculations  can  be  seen  to  significantly  over¬ 
estimate  the  ratio  aIR/aVIS  for  windspeeds  <  10  m/s. 

2. 2. 1.5  Naval  Research  Laboratory  (NRL)  Long-Path 
Infrared  Transmission  Measurements 

Long-path  atmospheric  transmission  measurements  were 
performed  by  NRL  during  the  period  1974-1979  using  laser- 
transmissometer-calibrated  high-resolution  Fourier  transform 
spectroscopy  [28,29],  During  this  measurement  program 
overwater  experiments  were  conducted  at  Cape  Canaveral  Air 
Force  Station  (CCAFS) ,  Florida  and  also  at  SNI  using  the 
same  optical  path  as  the  PMTC-OSP  measurements.  A  compen¬ 
dium  report  containing  a  tabulation  of  data  acquired  at 
CCAFS  was  published  in  1977  [30]  and  an  analysis  of  these 
measurements  will  be  presented  in  a  forthcoming  report 
[23] .  A  report  describing  the  NRL  long-path  transmission 
measurements  at  SNI  in  1979  will  be  available  in  the  near 
future  [24].  Selected  examples  of  analysis  drawn  from 
these  reports  are  available  and  will  be  considered  since 
they  contain  relevant  maritime  aerosol  infrared  extinction 
data. 
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FIGURE  16.  AEROSOL  IR  TO  VISIBLE  EXTINCTION  RATIOS  VERSUS 
12  HOURLY  AVERAGE  WINDSPEED.  A  COMPARISON  BE¬ 
TWEEN  SNI  EXPERIMENTAL  DATA,  NPGS  NORTH  ATLANTIC 
AEROSOL  DISTRIBUTION  DATA  AND  SUBSEQUENT  EXTINC¬ 
TION  RATIO  CALCULATIONS,  AND  THE  K/R  AND  LOWTRAN 
3B  MARITIME  MODEL  RESULTS. 
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Table  1  which  is  based  on  reference  23  contains  a 
summary  of  the  daily  observations  of  wind  speed  and  direc¬ 
tion,  visibility,  absolute  humidity,  air  temperature  and 
apparent  aerosol  extinction  (AAE)  values  determined  for  the 
different  laser  wavelengths  measured  in  the  1977  NRL-CCAFS 
experiment.  The  AAE  values  correspond  to  measured  extinc¬ 
tion  coefficients  for  HeNe  and  Nd-YAG  laser  wavelengths  and 
to  the  quantity  EXT-CMA  for  DF  and  C02  laser  wavelengths, 
where  EXT  are  measured  total  extinction  coefficients  and 
CMA  are  calculated  molecular  absorption  values  corresponding 
to  the  atmospheric  conditions  and  laser  wavelengths  studied 
during  a  particular  set  of  measurements.  When  several  laser 
frequencies  were  measured  during  a  particular  experimental 
run  (on  a  particular  day)  an  average  value  of  AAE  for  the 
several  measurements  was  determined.  Occasionally  long- 
path  extinction  measurements  for  many  laser  lines  were 
collected  (as  many  as  80  C02  laser  lines  on  some  days  and 
repeated  measurements  of  a  manifold  of  several  DF  laser 
lines) .  A  representative  value  for  AAE  appropriate  to  the 
spectral  interval  covered  by  the  several  laser  lines  was 
then  determined  and  used  as  the  value  appearing  in  Table  1. 
The  measurements  performed  for  overwater  wind  conditions 
(designated  by  OW  in  column  3  of  the  table)  are  considered 
to  be  most  representative  of  maritime  conditions. 

Ratios  of  selected  AAE  values  appearing  in  Table  1 
were  formed  and  are  compared  to  K/R  model  predictions  in 
Table  2.  The  ratio  of  the  experimental  data  for  the  dif¬ 
ferent  laser  wavelengths  are  designated  by  the  symbol  Re 
and  the  laser  identification  in  appropriate  columns  of  the 
table.  Data  were  selected  for  overwater  wind  conditions 
of  at  least  2  m/s  and  for  relative  humidities  >  65%.  One 
can  see  from  an  examination  of  the  entries  in  the  table 
that  the  K/R  model  predictions  for  Re  are  generally  lower 
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Dally  Record  of  Wind  Character,  Meteorological 
Parameters  and  Apparent  Aerosol  Extincnoc  (AAE)  Vali 
from  NRL  1977  CCAFS  Data  Base 
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Comparison  of  °ir/°vis  Values  from  NRL  1977 
CCAFS  Data  Base  to  Aerosol  Model  Predictions 
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than  those  observed  for  the  Nd-YAG  laser  transmission  data. 

On  the  averaae,  the  K/R  predictions  are  about  74%  of  the 
measured  values  for  Re  =  ajR/°VIS  for  the  Nc^-YAG  laser  data 
and  about  47%  of  the  values  measured  for  the  DF  laser  data. 
The  windspeed  values  correspondinq  to  these  data  are  general¬ 
ly  in  the  3  to  5  m/s  range  and  under  these  conditions,  as 
shown  in  Figure  2,  the  K/R  model  predicts  significantly 
lower  values  for  ajR/aVIS  than  does  either  version  of  LOW- 
TRAN. 


The  variability  in  the  data  for  both  wavelengths  does 
not  appear  to  be  particularly  well  represented  by  the  K/R 
model.  The  LOWTRAN  3B  and  LOWTRAN  5  models  show  better 
agreement  in  magnitude  of  Re  with  the  experimental  data 
than  does  the  K/R  model  but  the  scatter  in  Re  values  in 
this  limited  set  of  comparisons  is  too  large  for  this  small 
sample  set  to  provide  any  conclusive  results. 

Figures  17-21  shows  plots  of  DF  laser  extinction  data 
collected  during  the  NRL-SNI  measurements.  Measured  ex¬ 
tinction  coefficients  ( ot)  minus  calculated  molecular  ab¬ 
sorption  values  (CMA)  are  plotted  versus  wavenumber  for  the 
set  of  7  DF  laser  lines  used  in  the  experiment.  The  quan¬ 
tity  a  -  CMA  is  then  a  direct  measure  of  aerosol  attenuation 
over  the  4.07  km  path  provided  the  CMA  values  are  correct. 
Based  on  extensive  comparisons  of  similar  measurements  to 
molecular  absorption  models  [31,32],  the  CMA  values  are 
known  quite  well.  The  largest  remaining  uncertainty  in  the 
molecular  absorption  values  concerns  the  magnitude  of  the 
weak  absorption  due  to  the  water  vapor  continuum  in  the  DF 
laser  region. 

The  consistent  trends  showing  a  decrease  of  the  a  - 
CMA  values  versus  wavenumber  in  Figures  17-21  are  probably 
due  to  minor  discrepancies  in  the  spectral  shape  of  the 
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FIGURE  17.  MEASURED  EXTINCTION  COEFFICIENTS  MINUS  CALCULATED 
MOLECULAR  ABSORPTION  COEFFICIENTS  VERSUS  WAVE- 
NUMBER  FOR  NRL  MEASUREMENTS  AT  SNI  OF  SEVERAL  DF 
LASER  LINES.  FROM  TOP  TO  BOTTOM:  DATA  FOR  5/1/79 
AND  5/2/79  RESPECTIVELY. 
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FIGURE  18.  MEASURED  EXTINCTION  COEFFICIENTS  MINUS  CALCULATED 
MOLECULAR  ABSORPTION  COEFFICIENTS  VERSUS  WAVE- 
NUMBER  FOR  NRL  MEASUREMENTS  AT  SNI  OF  SEVERAL  DF 
LASER  LINES.  FROM  TOP  TO  BOTTOM:  DATA  FOR  5/2/79 
AND  5/3/79  RESPECTIVELY. 
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MEASURED  EXTINCTION  COEFFICIENTS  MINUS  CALCULATED 
MOLECULAR  ABSORPTION  COEFFICIENTS  VERSUS  WAVE- 
NUMBER  FOR  NRL  MEASUREMENTS  AT  SNI  OF  SEVERAL  DF 
LASER  LINES.  FROM  TOP  TO  BOTTOM:  DATA  FOR  5/8/79 
AND  5/9/79  RESPECTIVELY. 
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IGURE  21.  MEASURED  EXTINCTION  COEFFICIENTS  MINUS  CALCULATED  MOLECULAR 
ABSORPTION  COEFFICIENTS  VERSUS  WAVENUMBER  FOR  NRL  MEASURE¬ 
MENTS  AT  SNI  OF  SEVERAL  DF  LASER  LINES.  TOP-0826  PpT : 
BOTTOM- 101 8  PDT. 


1^0  continuum  absorption  model  used  in  generating  the  CMA 
values.  Average  values  of  a  -  CMA  were  determined  from 
several  of  the  curves  plotted  in  Figures  17-21  and  are 
tabulated  in  column  7  of  Table  3.  The  ratios  a__/a,r__  for 
the  experimental  data  are  tabulated  in  coluim  9.  Comparable 
ratios  determined  from  the  K/R  and  LOWTRAN  5  models  are 
tabulated  in  columns  10  and  11  respectively.  A  comparison 
of  the  entries  in  the  last  three  columns  of  the  table  shows 
that  the  experimental  values  (aDF/avis^  are  ^ower  than  either 
model  prediction  for  the  high  visibility  conditions  en¬ 
countered  on  5-1-77.  The  K/R  model,  although  predicting  a 
somewhat  lower  magnitude  than  experimentally  observed,  re¬ 
produces  the  trends  as  a  function  of  windspeed  shown  by 
the  data  better  than  the  LOWTRAN  5  model.  The  advantage 
of  the  windspeed  dependence  incorporated  into  the  K/R  model 
can  be  seen  from  this  comparison. 

2.2.2  MARINE  AEROSOL  DISTRIBUTION  MEASUREMENTS 

Comparisons  of  infrared  aerosol  extinction  values  de¬ 
rived  from  aerosol  spectrometer  measurements  in  maritime 
conditions  have  been  recently  performed  using  aircraft 
measurements  by  the  Naval  Ocean  Systems  Center  (NOSC)  [4] 
and  in  conjunction  with  British  weather  ship  operations  in 
the  North  Atlantic  by  the  Naval  Research  Laboratory  (NRL) 

[33].  Certain  NRL  measurements  performed  at  CCAFS  [30] 
were  supported  by  shore-based  aerosol  spectrometer  measure¬ 
ments.  During  analysis  of  the  NRL-CCAFS  data  [23]  it  be¬ 
came  apparent  that  for  overwater  wind  directions  aerosol 
extinction  values  derived  from  aerosol  spectrometer  data 
would  significantly  overestimate  the  actual  aerosol  extinc¬ 
tion. 

Figure  22  is  a  plot  of  differences  between  apparent 
aerosol  extinction  (AAE)  values  from  the  77  CCAFS  experiment 


49 


AAE-CAE  (km 


ow 

OL 

OW 

OL 

HeNe  0 

• 

DF  O 

♦ 

Nd-YAG  Q 

■ 

C02  V 

▼ 

FIGURE  22.  PLOT  OF  APPARENT  AEROSOL  EXTINCTION 
(AAE)  COEFFICIENTS  DERIVED  FROM  NRL 
LONG- PATH  TRANSMISSION  MEASUREMENTS 
AT  CCAFS  MINUS  CALCULATED  AEROSOL 
EXTINCTION  COEFFICIENTS  FOR  EACH  DAY 
OF  THE  EXPERIMENT. 
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and  values  calculated  from  shore-based  aerosol  spectrometer 
measurements  (calculated  aerosol  extinctions  -  CAE) .  Data 
for  the  HeNe  and  Nd-YAG  wavelengths  are  shown  in  the  left- 
hand  portion  of  the  figure  and  data  for  DF  and  CC>2  laser 
wavelengths  are  shown  on  the  right.  The  solid  symbols  cor¬ 
respond  to  measurements  during  conditions  of  overland  wind 
directions  (blowing  from  land  toward  the  ocean)  while  the 
open  symbols  represent  overwater  wind  conditions.  The  fact 
that  most  of  the  overwater  wind  differences  are  negative 
indicates  that  the  CAE  values  are  unrealistically  larae  on 
account  of  surf  generated  particles.  In  these  cases  the 
conditions  sampled  by  the  aerosol  spectrometer  were  not  re¬ 
presentative  of  the  overwater  optical  path.  This  problem 
has  also  been  observed  in  more  recent  measurements  at  SNI 
when  extinction  coefficients  derived  from  shore-based 
aerosol  spectrometer  measurements  were  compared  to  the  PMTC- 
OSP  transmis some ter  data. 

2. 2. 2.1  Naval  Ocean  Systems  Center  (NOSC)  Aerosol 
Spectrometer  Measurements 

NOSC  has  developed  an  airborne  aerosol  spectrometer 
measurement  capability  which  has  been  used  in  several  loca¬ 
tions  during  recent  years  and  notably  at  SNI  during  recent 
experiments.  Comparisons  of  vertical  profiles  of  extinc¬ 
tion  coefficients  at  3.75  ym  to  predictions  of  WGM,  K/R 
and  LOWTRAN  3B  models  were  recently  published  [4],  Figures 
presented  in  reference  4  show  that  the  WGM  model  generally 
predicts  higher  extinctions  and  consequently  shows  better 
agreement  with  the  measured  data  than  does  the  K/R  model. 

The  LOWTRAN  3B  model  does  not  compare  favorably  with  the 
experimental  data  since  it  shows  much  less  resemblence  to 
the  vertical  profile  of  c  than  does  either  of  the  other 
two  models.  Direct  comparisons  of  measured  and  calculated 
aerosol  distributions  shown  in  reference  4  show  that  the 
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K/R  predictions  more  closely  resemble  the  measured  distri¬ 
butions  for  larger  values  of  RH  than  measured,  possibly  in¬ 
dicating  that  the  particle  growth  factor  in  the  K/R  model 
should  be  increased.  Other  comparisons  in  this  same  report 
support  earlier  observations  that  the  K/R  model  underesti¬ 
mates  particle  concentration  and  consequently  aerosol  ex¬ 
tinction  at  lower  windspeeds. 

2. 2. 2. 2  Naval  Research  Laboratory  (NRL)  Aerosol 
Spectrometer  Measurements 

The  NRL  open-ocean  aerosol  spectrometer  data  were 
collected  during  a  28-day  station-keeping  session  aboard 
the  British  weather  ship  Fitzroy  in  late  June  1978.  As 
discussed  in  reference  32,  measurements  were  limited  to 
fair-weather  conditions. 

The  conditions  occurring  during  most  of  the  measure¬ 
ment  periods  were  windspeeds  between  3  and  6  m/s  and  RH 
values  between  80%  and  100%.  The  measurement  results  are 
presented  as  frequency  of  occurance  plots  in  reference  32 
and  show  that  the  most  frequently  occurring  value  of  aVIS 
(0.55  yin)  calculated  from  the  measured  aerosol  distribution 
data  was  ^  0.03  km  1.  Corresponding  values  for  3.8  ym  and 
10.0  ym  were  about  0.008  km  ^  and  0.003  km  ^  respectively. 

A  preliminary  tabulation  of  the  data  contained  in  reference 
32  was  available  prior  to  publication  [34]  and  was  used  in 
an  effort  to  derive  ratios  of  aerosol  scattering  in  marine 
air  at  3.75  and  10  ym  as  part  of  the  analysis  contained  in 
reference  23. 

Figure  23  is  a  plot  of  the  ratios  aiR/aVis  ^er^ve<* 
from  the  NRL  marine  aerosol  spectrometer  data.  As  evident 
from  the  distributions  of  data  points  seen  in  the  figure, 
a  large  amount  of  scatter  is  present  in  the  data.  Least 
squares  fits  to  the  groups  of  data  points  were  performed 
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FIGURE  23.  RATIO  OF  INFRARED  TO  VISIBLE  AEROSOL 

EXTINCTION  COEFFICIENTS  VERSUS  VISIBLE 
AEROSOL  EXTINCTION  COEFFICIENTS  DERIVED 
FROM  NRL  NORTH  ATLANTIC  AEROSOL  SPEC¬ 
TROMETER  DATA.  0:  FITTED  EQUATION 
IS  al . 06/ aVIS=0 .574  +  0.281  (a0.55), 
CORR.  COEFF  .  =  0.714;  <>  :  FITTED  EQUATION  IS 
a3. 8/o0.55=0.234  +  0.208  (a0.55),  CORR. 
COEFF. =0. 420?  □  :  FITTED  EQUATION  IS 
olO. 6/00.55=  -0.121  +  0.360  (o0.55), 
CORR.  COEFF. =0.597. 


and  the  fitted  equations  are  given  in  the  figure  legend. 

The  poor  correlation  coefficients  for  the  curve  fits  reflect 
the  fact  that  the  data  exhibit  a  large  amount  of  scatter. 

2.3  SUMMARY  COMPARISON  OF  MODEL  AND  MEASURED  MARITIME 

EXTINCTION  VALUES 

A  detailed  comparison  of  the  features  of  the  marine 
aerosol  extinction  models  has  been  discussed  in  Section  2.2. 
In  summary,  the  K/R  model  predicts  much  lower  values  of 
aiR//aVIS  for  low  windsPeeds  than  do  either  the  LOWTRAN  3B 
or  LOWTRAN  5  models.  The  three  models  predict  comparable 
results  in  the  range  of  windspeeds  around  10  m/s  and  for 
RH  values  between  70%  and  90%.  The  relative  humidity  de¬ 
pendence  of  LOWTRAN  5  is  more  severe  than  the  K/R  model 
above  and  below  these  values. 

Comparisons  of  experimental  long-path  aerosol  extinc¬ 
tion  values  show  that  the  seemingly  unsophisticated  LOWTRAN 
3B  model  represents  measured  values  reasonably  well  (see 
Figure  5) .  The  TNO  Laboratory  data  which  were  collected 
under  moderate  visibility  conditions  were  seen  to  be  more 
useful  for  aerosol  model  comparison  purposes  than  the 
Australian  DSTO  measurements  performed  under  very  high 
visibility  conditions. 

All  of  the  examples  of  aj^/aVIS  values  obtained  from 
optical  transmission  data  that  were  shown  in  Section  2.2 
exhibit  large  amounts  of  scatter  which  are  probably  attri¬ 
butable  to  the  influence  of  different  air  masses  (variations 
in  the  relative  concentrations  of  continental  and  maritime 
air)  . 

The  OPAQUE  data  (Figures  11  and  12)  show  typically 
large  amounts  of  scatter  and  coarse  agreement  with  the  SAI 
water  haze  model  for  3-5  urn  but  not  for  8-12  urn  where  the 
ratio  aIR/aVIS  is  nearly  always  lower  than  the  LOWTRAN 
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model  predictions.  The  limited  examples  of  OSP  data  shown 
in  Figures  13-15  are  intriguing  in  that  they  indicate  that 
a  simple  power-law  dependence  of  a  on  a  agrees  with 
the  data.  Presumably  the  data  relate  to  a  nearly  constant 
air  mass  and  we  might  expect  to  see  different  slopes  on 
graphs  of  this  same  type  for  different  air  masses. 

The  importance  of  long  averaging  times  for  wind  speed 
values  used  in  modeling  aerosol  extinction  ratios  is  demon¬ 
strated  by  the  comparisons  shown  in  Figure  16  of  the  PMTC- 
SNI  data  to  the  K/R  model.  Reasonably  good  agreement  is 
seen  to  be  obtained  by  the  use  of  extended  time  average 
windspeeds . 

The  NRL  long-path  optical  data  contained  in  Tables  1- 
2  and  Figures  17-21  and  Table  3  show  no  obvious  agreement 
with  any  of  the  three  aerosol  models  considered  although 
the  K/R  model  shows  better  agreement  at  low  values  of  wind- 
speed  with  the  NRL-CCAFS  data  (Table  2)  than  does  either 
version  of  LOWTRAN. 

The  NRL-SNI  data  were  collected  during  conditions  that 
conform  more  nearly  to  the  region  of  greatest  overlap  in 
the  model  predictions,  i.e.,  ^  10  m/s  windspeeds  and  80% 

RH.  On  the  average  the  K/R  model  seems  to  predict  aIR/crVIS 

values  closer  to  observations  than  does  the  LOWTRAN  5  model 

The  aerosol  spectrometer  measurements  performed  by 
NOSC  [4]  indicate  that  the  K/R  model  RH  dependence  may  be 
too  weak  and, as  seen  in  other  comparisons,  that  better  agree 
ment  with  measured  values  of  the  ratio  of  infrared  to  vi¬ 
sible  scattering  occurs  for  windspeeds  around  10  m/s. 

The  NRL  weather  ship  aerosol  spectrometer  data  are 
seen  to  exhibit  large  amounts  of  scatter  in  the  ratios 
alR/aVIS'  emP^as^-zin9  the  ambiguity  inherent  in  comparisons 
of  infrared  aerosol  extinction  data  and  models. 
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2.4  ANTICIPATED  RANGES  OF  MARINE  ATMOSPHERIC  EXTINCTION 

AND  BACKSCATTER 

A  summary  of  the  ranges  of  aerosol  scattering  values 
for  the  various  measurements  described  in  Section  2.3  is 
listed  in  Table  4.  Corresponding  ratios  of  aT_/a.,Te  for 
1.06  ym,  3.8  ym  and  10.6  ym  derived  from  the  extremes  of 
these  data  are  also  included.  As  expected,  the  median  values 
of  the  ratio  aIR/aVIS  decrease  with  increasing  wavelength. 

The  very  low  visibility  OPAQUE  data  exhibit  quite  large 
values  for  a  and  g  which  are  not  seen  in  the  other 
data.  The  range  of  values  of  aIR/aVIS  for  1.06  ym  is  about 
a  factor  of  2,  increasing  to  a  factor  of  about  4  for  3.8  ym 
and  as  much  as  a  factor  of  7  for  the  SNI-K/R  data  at  10.6  ym. 

Such  an  increase  in  the  range  of  values  of  a  /a  is 

IR  VIS 

in  keeping  the  increased  variability  of  particle  density 
for  large  sizes  relative  to  the  concentrations  of  smaller 
particles  effective  in  scattering  shorter  wavelengths. 

We  have  discussed  the  variations  of  model-predicted 
ratios  of  infrared  to  visible  aerosol  extinction  and  compared 
measured  values  to  model  predictions  in  Sections  2.2  and  2.3 
and  found  no  obvious  "best  overall  representation"  by  any 
one  of  the  models.  It  thus  seems  reasonable  to  exercise  one 
or  two  of  the  models  for  a  typical  set  of  maritime  condi¬ 
tions  for  the  purposes  of  estimating  a  LIDAR  system  perfor¬ 
mance  at  sea,  and  then  to  evaluate  the  excursions  from  this 
"nominal  case"  based  on  the  model  and  data  review  presented 
in  the  previous  sections. 

Table  5  contains  the  results  of  calculations  performed 
with  the  maritime  component  of  the  LOWTRAN  5  model.  Aerosol 
extinction,  absorption,  scattering  and  backscattering  were 
calculated  for  the  four  laser  wavelengths  shown  and  for  2  km 
and  5  km  visibilities  and  for  70%  and  95%  RH. 
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Summary  of  Observed  Aerosol  Extinction 
Values  and  Ratios  of  oTn/o,.TO 

IK  VIb 
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Calculations  of  Aerosol  Extinction,  Absorption,  Scattering 
and  Backscatter  for  Different  Wavelengths, 
Visibilities  and  Relative  Humidities  Using  The 
LOWTRAN  5  Maritime  Model  (Oceanic  Component) 
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We  are  interested  in  how  the  ratio  of  scattering  to 
extinction  behaves  as  a  function  of  wavelength  since  it  is 
required  that  this  ratio  remain  constant  if  we  wish  to 
implement  the  "Klett  method"  [35,36]  for  processing  LIDAR 
return  signals.  The  Klett  method  will  be  discussed  in 
greater  detail  in  Section  3.2. 

One  sees  that  the  ratio  of  scattering  to  extinction 
does  not  change  drastically  between  0.55  ym  and  3.8  um 
wavelengths  but  that  it  is  appreciably  smaller  at  10.6  um 
due  to  increased  absorption  by  the  water  droplet  aerosols 
at  the  longer  wavelength.  The  70%  RH  cases  for  2  km  and 
5  km  visibilities  are  seen  to  be  somewhat  worse  than  the 
95%  RH  case.  The  ratio  of  scattering  to  extinction  changes 
by  about  20%  between  0.55  um  and  3.8  um  in  the  former  case 
and  only  by  about  6%  in  the  latter  case. 

The  calculated  backscatter  ratio  (km  1  sterad  1)  list¬ 
ed  in  the  sixth  column  of  Table  5  is  required  to  estimate 
the  LIDAR  return  signal  strength.  One  can  see  that  both 
the  scattering  magnitude  a  and  the  backscatter  coefficient 
8  decrease  with  increasing  wavelength. 

The  backscatter  coefficient  6  is  determined  by  per¬ 
forming  a  Mie  angular  scattering  calculation  based  on  the 
marine  aerosol  distribution  generated  by  the  particular 
model  being  used,  e.g.  LOWTRAN,  K/R  or  WGM.  The  product 
of  the  angular  scattering  phase  function  for  a  scattering 
angle  of  180°  and  the  magnitude  of  the  scattering  coeffi¬ 
cient  then  determines  the  value  of  6  for  particular  scat¬ 
tering  condition.  Examples  of  the  scattering  phase  function 
for  an  oceanic  aerosol  for  1.06  um  and  10.6  um  are  shown  in 
Figures  24  and  25.  One  can  see  that  scattering  at  1.06  um 
for  this  type  of  aerosol  is  sharply  peaked  in  the  forward 
direction  while  it  is  more  nearly  anqularly  isotropic  at 
10.6  ym.  Even  though  the  magnitude  of  o  is  much  lower  at 


10.6  ym  than  at  1.06  ym  for  the  oceanic  aerosol  the  ratio 
of  S/ct  is  actually  larger  at  the  longer  wavelength. 

Similar  results  to  those  tabulated  in  Table  5  and 
shown  in  Figures  24  and  25  were  performed  for  the  continen¬ 
tal  component  in  LOWTRAN  5.  The  angular  scattering  phase 
function  results  at  10.6  ym  are  quite  similar  to  those 
shown  in  Figure  25  but  the  results  for  the  continental 
component  at  10.6  ym  do  not  show  the  strong  forward  scat¬ 
tering  peak. 

Figure  26  is  a  plot  of  the  LOWTRAN  5  continental  com¬ 
ponent  angular  scattering  phase  function  for  1.06  ym  show¬ 
ing  a  more  uniform,  nearly  angular  isotropic  behavior. 

Calculations  were  also  performed  using  the  WGM  model 
for  70%  and  95%  RH,  for  windspeeds  between  0  and  10  m/s 
and  for  0.55  ym,  3.8  ym  and  10.6  ym  wavelengths.  The  re¬ 
sults  of  these  calculations  are  shown  in  Table  6.  Here  we 
see  that  the  ratio  of  scattering  to  extinction  is  even  more 
constant  between  0.55  ym  and  3.8  ym  wavelengths  than  for 
the  LOWTRAN  5  oceanic  component  case.  However  we  see  that 
this  ratio  falls  from  a  value  near  unity  for  3.8  ym  to  about 
0.4  to  0.6  at  10.6  ym,  raising  doubts  about  the  validity  of 
the  assumptions  required  by  the  Klett  inversion  method  at 
10.6  ym. 

Results  of  calculations  of  the  backscatter  coefficient 
8  for  the  oceanic  and  rural  components  of  LOWTRAN  5  for  the 
two  cases:  a)  2  km  visibility  and  95%  RH  and  b)  5  km  visi¬ 
bility  and  70%  RH  are  plotted  as  a  function  of  wavelength 
in  Figure  27.  Results  of  calculations  using  an  advection 
fog  model  for  2  km  and  5  km  visibilities  [22]  are  also 
shown  in  the  figure. 

When  the  results  such  as  those  shown  in  Table  5  and 
Figures  24-26  for  continental  and  oceanic  aerosol  components 
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FIGURE  24.  OCEANIC  AEROSOL  SCATTERING  PHASE  FUNCTION  FOR  1.06  ym 
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FIGURE  25.  OCEANIC  AEROSOL  SCATTERING  PHASE  FUNCTION  FOR  10.59  pm 
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FIGURE  26.  CONTINENTAL  AEROSOL  SCATTERING  PHASE  FUNCTION  FOR  1.06  ym. 
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FIGURE  27.  CALCULATED  AEROSOL  BACKSCATTER  COEFFICIENTS 

VERSUS  WAVELENGTH  FOR  DIFFERENT  AEROSOL  MODELS 
WITH  DIFFERENT  VISIBILITIES  AND  RELATIVE 
HUMIDITIES . 
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are  combined  compensating  effects  occur.  An  actual  marine 
aerosol  distribution  will  consist  of  certain  proportions 
of  one  component  to  the  other.  The  behavior  of  quantities 
like  the  scattering  phase  function  (e.g.  Figures  24  and  26) 
for  different  components  when  combined  in  a  given  ratio 
will  produce  results  which  are  dependent  upon  the  relative 
proportions  of  continental  and  marine  components.  That  is, 
results  of  calculations  of  the  ratio  of  6/a  as  a  function 
of  wavelength  will  depend  upon  the  relative  proportions  of 
the  two  components  chosen.  Figure  28  is  an  example  of  a 
plot  of  the  ratio  6/a  versus  wavelength  for  LOWTRAN  5  using 
equal  components  of  rural  and  oceanic  aerosol  distributions 
The  5  km  visibility  case  is  shown.  Results  for  the  5  km  vi 
sibility  fog  model  are  shown  also. 

Using  the  results  shown  in  Figure  28  we  can  perform 
estimates  of  backscattered  signal  levels  given  a  particular 
LIDAR  system  configuration.  The  results  of  these  calcula¬ 
tions  will  represent  a  "nominal  case"  for  operation  in  a 
maritime  atmosphere.  The  information  discussed  in  Sections 
2.2  and  2.3  can  then  be  used  as  a  basis  for  interpreting 
the  nominal  case  results  and  estimating  the  changes  in  sys¬ 
tem  performance  that  might  result  from  excursions  in  scat¬ 
tering  and  backscatter  corresponding  to  the  range  of  in¬ 
frared  and  visible  scattering  variations  documented  earlier 

A  brief  review  of  LIDAR  measurement  techniques  and 
approaches  will  be  presented  in  Section  3  and  some  of  the 
more  important  considerations  involved  in  implementing  a 
LIDAR  measurement  system  on  shipboard  *fill  be  discussed  in 
Section  4. 


68 


i-ava3is)  -% 


M  yum) 


FIGURE  28.  CALCULATED  RATIO  OF  AEROSOL  BACKSCATTER  TO 
EXTINCTION  COEFFICIENTS  VERSUS  WAVELENGTH 
FOR  DIFFERENT  AEROSOL  MODELS  WITH  DIFFERENT 
VISIBILITIES  AND  RELATIVE  HUMIDITIES. 
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LIDAR  MEASUREMENT  SYSTEM  APPROACHES 

3.1  TECHNIQUES  USING  CW  LASERS 

In  an  attempt  to  circumvent  the  problems  associated 
with  conventional  pulsed  LIDAR  systems,  several  variations 
using  CW  lasers  have  been  proposed  and  investigated  [37-43]. 
Reduction  of  eye  safety  problems  realized  by  the  use  of  low 
power  CW  laser  sources  as  opposed  to  high  peak  power  pulsed 
lasers  is  an  attractive  feature  of  CW  LIDARS.  The  costs 
and  complexity  associated  with  high-speed  digital  process¬ 
ing  hardware  can  be  reduced  by  using  simplified  analog  data 
acquisition  and  processing  of  CW  signals.  The  CW  LIDAR 
techniques  which  have  been  investigated  and  are  described 
in  references  37-43  can  be  divided  into  three  types  of 
systems,  namely  a)  those  using  a  fixed  frequency  of  modu¬ 
lation  of  the  CW  laser  source  (chopped  CW  source) ,  b)  sys¬ 
tems  where  the  CW  laser  modulation  frequency  is  swept  thru 
a  range  ("chirped"  modulation  of  the  CW  source,  analogous 
to  a  chirped  radar  system)  and  c)  systems  where  range  dis¬ 
crimination  is  accomplished  by  spatial  analysis  of  an  image 
formed  by  the  backscattered  return  light.  The  beam- image- 
profile  approach  described  in  references  40-42  is  an  example 
of  this  type  of  system.  Considerations  involved  in  these 
types  of  CW  LIDAR  systems  will  now  be  briefly  reviewed. 

3.1.1  FIXED  FREQUENCY  MODULATED  LIDARS  (CHOPPED  CW 
LASER  SOURCES) 

Stokes  has  reviewed  the  work  of  Kreid  [37],  Paulson 
[38],  and  Bufton  and  Iyer  [39]  in  a  recent  report  [40], 

Kreid  presented  an  engineering  evaluation  but  did  not  ex¬ 
perimentally  verify  his  results.  The  other  authors  reported 
on  the  results  of  measurements  using  a  CW  LIDAR.  The  basic 
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approach  in  these  systems  is  a  comparison  of  the  phase  of 
the  backscattered  return  signal  to  the  phase  of  the  source. 
This  can  be  seen  from  the  following  analysis. 

The  LIDAR  return  at  time  t  given  by  the  LIDAR  equation 
is : 

P(t)  =  /* 3(2)  Ad  exp(-/”2a(C)dC)d2  (1) 

z 

where  8(z)  is  the  backscatter  coefficient,  a  is  the  extinc¬ 
tion  coefficient,  z  is  range,  PQ  is  the  modulated  laser 
power,  c  is  the  velocity  of  light,  PQ(t)  =  PQ  +  sin(wt), 
Ad  is  the  collection  optic  area  and  f(z)  is  the  correction 
for  actual  receiver  collection  cone.  For  an  exact  result, 
the  atmosphere  must  be  homogeneous  so  that: 

a  *  a  (z)  ,  8  =  6  (z)  .  (2) 


If  we  also  substitute  the  expression  for  PQ(t)  we  obtain: 


P(t)  *  Ad6  p  exp  (-2az)  f  [DC  term]  + 


AdB  %Qr0  exp  -2az  sinw  (t-2z/c)  f -  ■  dz  [AC  term] 


If  we  set  the  AC  term  equal  to: 


(3) 


P(t)  »  A  sinwt  -  B  coswt 

=  (A2  +  B2)1*  sin(wt  -  <f>)  , 


(4) 
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then  A  &  B  may  be  evaluated: 


l\j  oo 

A  =  V  Vo 

exp (-2az) 

cos  (2ajz/c) 

f (z) dz 

2 

z 

(5) 

b  -  v  V" 

exp (-2az ) 

sin  (2coz/c) 

f  (z )  dz 

2 

(6) 

The  phase  angle  4>  is  then  evaluated  as : 

<j>(a,w)  =  tan-1  (B/A)  (7) 

and  from  this  the  extinction  coefficient  can  be  determined. 

A  CW  LIDAR  system  was  constructed  and  tested  by  Paulson 
who  pointed  out  that  a  visibility  measurement  with  such  a 
system  is  most  sensitive  to  atmospheric  conditions  near  the 
convergence  distance  (point  of  crossover  of  the  LIDAR  trans¬ 
mitted  beam  and  receiver  f ield-of-view) . 

In  order  to  investigate  the  CW  LIDAR  phase  shift  sen¬ 
sitivity  as  a  function  of  crossover  distance  and  atmospheric 
homogenity,  a  simulation  of  a  CW  LIDAR  response  was  per¬ 
formed  for  several  cases.  Results  of  this  simulation  are 
shown  in  Figure  29  where  the  phase  angle  response  of  the 
system  is  plotted  against  visibility. 

The  response  of  the  system  to  a  uniform  atmosphere  of 
the  indicated  visibility  is  shown  by  the  solid  curves  for 
values  of  the  LIDAR  crossover  range,  Zo  =  10  m,  100  m  and 
200  m.  These  curves  may  be  considered  as  calibration  curves 

Response  of  the  LIDAR  to  a  non-uniform  atmosphere  is 
shown  by  the  dashed  curves  for  the  same  values  of  Zo.  Here 
the  LIDAR  is  located  in  a  clear  air  region  and  a  region  of 
uniform  aerosol  is  located  a  distance  Z'  away  from  the  LIDAR 
For  the  case  of  Zo  >>  Z'  where  the  system  crossover  occurs 
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FIGURE  29.  CALCULATED  PHASE  ANGLE  OF  CV7  LIDAR  RESPONSE  VERSUS 
VISIBILITY  FOR  UNIFORM  AND  NON-UNIFORM  ATMOSPHERIC 
AEROSOL  DISTRIBUTIONS. 
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well  beyond  the  atmospheric  discontinuity,  the  uniform  at¬ 
mosphere  calibration  response  curve  and  the  inhomogeneous 
atmosphere  response  curve  essentially  coincide  showing  that 
the  system  performs  quite  well. 

The  response  curve (s)  also  shows  that  for  Zo  =  200  m 
the  sensitivity  of  phase  angle  response  to  changes  in  visi¬ 
bility  is  good  and  is  only  starting  to  flatten  out  for  R^ 
j>  10  km,  indicating  that  useful  measurements  can  be  obtained 
for  this  and  higher  visibilities.  For  cases  where  the  cross¬ 
over  range  lies  near  the  boundary  of  the  atmospheric  inhomo- 
genity  (Zo  =  Z')  or  inside  it  (Zo  <<  Z')  the  system  response 
to  the  inhomogoneous  atmosphere  differs  from  that  in  a  uniform 
atmosphere.  The  middle  and  lower  curves  in  Figure  29  de¬ 
monstrate  these  cases.  With  Zo  and  Z'  each  set  equal  to 
100  m  we  see  that  reasonable  results  are  still  obtained 
although  the  visibility'  is  slightly  underestimated  in  a 
non-uniform  atmosphere.  The  lower  curves  show  that  visi¬ 
bility  estimates  can  be  in  error  by  20%  to  50 a  for  the  case 
where  for  the  Zo  <<  Z'.  The  use  of  return  beam  phase  com¬ 
parisons  as  a  measure  of  visibility  with  a  CW  LIDAR  system 
is  seen  to  be  sensitive  to  the  relative  locations  of  the 
system  optical  crossover  and  atmospheric  inhomogeneities. 

The  location  of  the  optical  crossover  with  such  a  system  is 
very  sensitive  to  system  optical  alignment.  This  fact  can 
be  inferred  from  Figure  29  by  noting  that  the  three  "cali¬ 
bration-response"  curves  for  a  uniform  atmosphere  differ 
appreciably  depending  upon  the  location  of  the  optical  cross¬ 
over.  Paulson  [38]  and  Bufton  and  Iyer  [39]  have  also  re¬ 
ported  that  optical  alignment  sensitivity  is  the  principal 
source  of  error  in  the  operation  of  such  systems. 
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3.1.2  SWEPT  FREQUENCY  MODULATED  LIDARS  ("CHIRPED"  CW 
LASER  SOURCES) 

A  LIDAR  based  on  f requency-swept-modulation  of  the 
source,  analogous  to  a  "chirped"  radar,  is  described  by 
Stokes  [40] .  The  analysis  and  evaluation  of  this  approach 
presented  by  Stokes  shows  that  the  system  is  primarily  sen¬ 
sitive  to  nonuniformities  and  discontinuities  in  the  atmos¬ 
phere  and  is  not  a  practical  approach  for  measuring  extinc¬ 
tion  due  to  uniformly  distributed  scatters.  The  FM  techni¬ 
que  has  been  demonstrated  to  provide  range  resolved  infor¬ 
mation  about  inhomogeneities,  but  it  does  not  respond  to 
the  average  backscatter  and  hence  is  not  useful  for  measuring 
visibility. 

It  is  interesting  to  speculate  on  the  possibility  that 
an  FM-CW  LIDAR  could  also  be  operated  as  a  fixed- frequency 
modulated  system  using  phase  discrimination.  The  comple¬ 
mentary  character  of  the  information  provided  by  the  two 
types  of  systems  would  then  provide  all  of  the  required 
information.  The  FM-CW  LIDAR  return  could  be  used  to  pro¬ 
vide  the  map  of  inhomogeneities  required  to  properly  inter¬ 
pret  the  CW  LIDAR  return. 

3.1.3  BEAM  IMAGE  PROFILING  (BIP)  LIDAR 

The  BIP  LIDAR  technique  was  developed  at  NOSC  and  is 
also  described  by  Stokes  [40].  This  is  a  simplified  CW 
LIDAR  approach  that  uses  geometry  rather  than  time  to  se¬ 
parate  the  signal  return  from  different  ranges.  The  tech¬ 
nique  has  been  studied  and  under  development  for  some  time 
at  NOSC  [40-43].  The  offset  in  the  optical  axes  of  the  BIP 
LIDAR  laser  transmitter  and  receiver  causes  the  backscatter- 
ed  return  signal  to  form  a  line  image  in  the  receiver  tele¬ 
scope  at  an  oblique  angle  to  the  receiver  telscope  axis. 

In  the  BIP  LIDAR  the  distance  of  a  point  in  the  line  image 
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from  the  telescope  axis  is  inversely  proportional  range. 

The  mapping  of  the  light  backscattered  from  the  transmitted 

beam  along  the  length  of  the  image  goes  inversely  as  the 

square  of  the  range,  building  in  an  automatic  compensation 
2 

for  the  1/r  fall-off  of  a  LIDAR  return  signal.  Several 
discreet  detectors  and  signal  processing  channels  can  be 
used  to  measure  the  image  brightness  at  different  locations 
along  the  image  and  thus  measure  the  return  signal  from 
different  ranges.  Recent  developments  in  commercially 
'available  linear  Si  diode  array  detectors  containing  as 
many  as  1024  detectors  and  packaged  in  integrated  circuit 
housings  along  with  integral  clocking  electronics  provide 
a  detector  module  ideally  suited  to  this  application. 

Jensen  [43]  has  summarized  the  results  of  tests  and  inter- 
compa  'isons  of  a  BIP  LIDAR  system  with  simultaneous  measure¬ 
ments  using  other  instruments  at  NOSC.  The  conclusions 
reached  were  that  intercomparisons  of  the  BIP  LIDAR  measure¬ 
ments  with  telephotometers,  nephelometers  and  aerosol  spec¬ 
trometers  generally  showed  agreement  to  within  a  factor  of 
two.  The  critical  nature  of  system  optical  alignment  and 
the  dependence  of  the  system  calibration  upon  the  thermal 
stability  of  this  alignment  seem  to  be  the  major  factors 
limiting  the  usefulness  of  the  BIP  system  for  use  in  routine 
visibility  measurements.  A  three-dimensional  computer 
search  routine  was  required  to  determine  the  image  vanish¬ 
ing  point  on  the  detector  array  before  data  could  be  cali¬ 
brated  in  visibility. 

3.1.4  SUMMARY  OF  ADVANTAGES  AND  DISADVANTAGES  OF  CW 
LIDAR  APPROACHES 

There  are  several  potential  advantages  to  the  use  of 
CW  LIDAR  techniques  in  designing  a  single-ended  visibility 
measurement  system  which  make  the  evaluation  of  such  systems 
for  this  application  worthwhile.  The  same  arguments  pro 
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and  con  generally  apply  to  the  three  variations  of  CW  systems 
previously  described.  Some  important  advantages  of  a  CW  system  are 

•  Eye  safety  problems  are  simplified  due  to  the  use 
of  a  CW  laser. 

•  Low  electro-magnetic-interference  (EMI)  generation. 

•  The  assumed  relationship  between  extinction  and 
backscatter  is  not  required  (as  in  the  application 
of  the  Klett  method  with  a  pulsed  LIDAR) . 

•  Some  simplification  in  electronic  signal  processing 
complexity  is  possible  since  high-speed  analog- 
to-digital  converters  are  not  required.  The  signal 
processing  (e.g.  phase  comparison  of  transmitted 
and  backscattered  beams)  can  be  readily  implemented 
with  analog  electronics. 

2 

•  Inherent  optical-geometrical  compensation  for  1/r 
signal  loss  (BIP  LIDAR  only) . 

The  findings  presented  in  references  37-43  together 
with  the  analysis  shown  in  Figure  29  suggest  several  limi¬ 
tations  inherent  in  the  CW  LIDAR  approach  when  these  systems 
are  considered  for  routine  use  in  a  visibility  measurement 
application.  Principal  among  these  disadvantages  are  the 
following: 


•  Difficult  optical  alignment  requirements  limit  the 
utility  and  measurement  accuracy  with  all  of  the 
systems  described  in  the  literature.  This  is  the 
overriding  limitation  of  these  systems. 


•  The  phase  comparison  method  using  a  fixed-frequency- 
modulated  CW  source  can  only  work  exactly  for  a 
uniform  atmospheric  path. 

•  The  swept-frequency-modulated  "chirped"  system  has 
been  shown  to  respond  only  to  discontinuities  in 
backscatter  and  hence  is  not  useful  for  visibility 
measurements  in  a  homogeneous  atmosphere. 
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•  The  BIP  LIDAR  system  has  an  inherent  non-linear 
range  resolution  due  to  the  compression  of  the  re¬ 
turn  signal  from  distant  ranges  near  the  image 
vanishing  point. 

3.2  APPLICATION  OF  PULSED  LIDAR  SYSTEMS  TO  VISIBILITY 

MEASUREMENTS 

3.2.1  GENERAL  CHARACTERISTICS  OF  A  PULSED  LIDAR 
SYSTEM 

A  LIDAR  using  a  pulsed  laser  source  is  generally  much 
simpler  in  optical  design  and  construction  than  the  CW  LIDAR 
systems  discussed  earlier  and  described  in  references  36- 
42.  For  shipboard  visibility  (aerosol  extinction)  measure¬ 
ment  applications  a  monostatic  LIDAR  system  would  be  most 
practical.  In  this  configuration  the  laser  source  and  re¬ 
ceiver  are  nearly  co-located  adjacent  to  one  another  or 
possibly  coaxial  and  the  projected  beam  and  receiver  field- 
of-view  converge  or  overlap  at  a  distance  of  several  tens 

of  meters  away  from  the  LIDAR.  A  relatively  short  laser 

_o 

pulse  of  about  10  sec  duration  is  used  in  order  to  obtain 
good  range  resolution  of  the  LIDAR  backscatter  profile. 

Some  means  is  required  to  record  and  store  the  rapidly 
varying  time  dependent  backscatter  signal  so  that  the  sig¬ 
nal  strength  as  a  function  of  time,  i.e.,  range  may  be 
determined.  A  transient  recording  capability  or  very  high 
speed  analog- to-digital  converter  is  required  to  capture 
the  high-speed  transient  LIDAR  return  signal.  These  con¬ 
siderations  will  be  discussed  in  greater  detail  in  Section 
3.3. 

3.2.2  SIGNAL  PROCESSING  METHODS 

The  LIDAR  equation  for  the  instantaneous  received 
power  P(r)  can  be  written: 
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P  (r ) 


(8) 


-  Po  ¥  AR  ^-T1  exPl'2^  o  (r '  )dr '  1 
r 

where  PQ  is  the  transmitted  power  at  time  tQ,  c  is  the 
velocity  of  light,  t  is  the  pulse  duration,  AD  is  the  ef- 
fective  receiver  area,  r  is  the  range,  B(r)  is  the  volume 
backscatter  coefficient  and  a(r)  is  the  attenuation  coef¬ 
ficient.  The  LIDAR  equation  can  be  re-written  in  terms  of 
2 

the  1/r  compensated,  logarithmic  signal  variable  S(r) : 

S(r)  =  In (r2P (r) ) .  (9) 

The  LIDAR  equation  (8)  then  becomes: 

S(r)  -  S (r  )  =  In  8/B.  -  2 fr  adr'  (10) 

ro 

where  B0  is  evaluated  at  the  reference  range  rQ. 

3. 2. 2.1  Slope  and  Ratio  Methods  of  Signal  Inversion 

The  "slope"  method  requires  that  the  atmosphere  be 
homogeneous  so  that : 


dB  _ 

a?  ~ 


0 


In  this  case,  the  volume  extinction  is  simply: 


a 


1  dS 

=  _  7  3F 


(11) 


(12) 


Since  in  reality  S(r)  is  never  exactly  a  linear  function 

of  r,  i.e.,  a  is  never  exactly  independent  of  range,  a 

least-squares  fitting  procedure  may  be  applied  to  evaluate 
dS 

gp.  In  doing  so,  one  hopes  that  path  discontinuities  will 


80 


simply  average  out.  A  variant  of  the  "slope"  method  is  to 
break  up  the  return  signal  into  short  segments.  In  this 
way,  the  probability  of  encountering  a  scattering  discon¬ 
tinuity  within  a  given  segment  is  decreased  and  the  as¬ 
sociated  probability  that  ^  is  small  or  zero  is  increased. 
When  summing  over  a  path,  it  is  hoped  that  one  can  identify 
discontinuities  and  correct  them.  This  approach  works  best 
when  changes  in  scattering  occur  at  sharp  boundaries  so  that 
their  effect  can  be  detected  and  ignored. 

The  "ratio"  method  is  essentially  identical  to  the 
previous  method  in  assumptions  and  results.  Its  formula¬ 
tion,  however,  points  out  the  fact  that  if  two  points  in 
a  LIDAR  return  share  the  same  backscatter  coefficients, 
the  transmission  along  the  path  joining  those  points  can 
be  evaluated  independent  of  the  scattering  behavior  along 
the  path.  Thus,  the  ratio  technique  is  useful  for  evalu¬ 
ating  the  transmission  between  two  points  which  bound  the 
region  of  inhomogeneity  in  an  otherwise  homogeneous  atmos¬ 
phere  . 


3.2. 2. 2  Klett  Method  of  Signal  Inversion 

The  most  general  inversion  technique  is  a  solution  to 
the  LIDAR  equation  assuming  a  power  law  relationship  be¬ 
tween  volume  extinction  and  backscatter: 

6  =  •  (13) 


Then  the  LIDAR  equation  becomes: 


dS  _  k  do 
3r  ~  o  3r 


2a 


whose  solution  can  be  written  as : 


(14) 
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S  (r) 


S  (r) 


) dr  1 .  (15) 


air)  =  exp 


2 


,r 

fr  exp 
o 


-  S(ro} 


This  solution  allows  arbitrary  path  structure  so  long  as  the 
power  law  relation  between  a  and  B  is  valid.  Generally, 
this  equation  has  not  been  used  due  to  its  extreme  sensiti¬ 
vity  to  noise  and  the  initial  determination  of  oQ  at  range 
rQ  (evaluated  using  the  slope  method).  Klett  [35,36]  noted 
in  his  paper  that  the  prime  reason  for  this  effect  is  the 
mathematical  instability  in  the  equation  structure.  Klett 
shows  that  the  reason  for  this  is  that  cq  is  typically 
evaluated  near  the  LIDAR  itself  so  that  a(r)  is  always 
evaluated  at  ranges  where  the  return  signal  and  thus  the 
S/N  is  diminished.  It  is  the  differencing  of  two  rather 
large  numbers  to  yield  a  small  result  in  the  denominator 
that  cause  the  trouble.  This  situation  is  greatly  reduced 
if  aQ  is  evaluated  at  a  range  rm  greater  than  those  of  in¬ 
terest.  Then  the  LIDAR  equation  becomes: 


a(r)  =  exp(- 


S  (r)  -  S  (r  ) 


m 


-yv + 1  flm  exp< 


Sir)  -  S  (r  ) 
m 


)  dr ' 


(16) 


and  stability  of  the  solution  is  greatly  improved  by  this 
simple  modification. 

3.3  CHARACTERISTICS  OF  THE  ANGVS/5  MODEL-1  VISIOCEILOMETER 

The  Army  has  been  quite  successful  in  developing  a 
hand-held  LIDAR  which  is  suited  to  Army  requirements.  The 
development  was  initiated  at  the  Night  Vision  and  Electro- 
Optics  Laboratory  (NVEOL)  [44]  and  is  an  on-going  activity  at 
the  Atmospheric  Sciences  Laboratory  (ASL)  at  the  White 
Sands  Missile  Range  (WSMR)  [45] . 
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The  LIDAR  is  a  visioceilometer  (ASL  Model-1)  based  on 
the  AN/GVS-5  portable  laser  rangefinder.  This  device  is 
quite  small  —  occupying  a  volume  less  than  one  cubic  foot. 
A  second  generation  (Model-2)  version  will  be  designed  for 
hand-held  operation  and  will  be  battery  powered  using  a 
belt  pack.  Each  module  of  the  Model-2  version  is  expected 
to  weigh  between  5  and  6  pounds .  The  ASL  development  pro¬ 
gram  has  strongly  emphasized  miniaturization  since  the 
design  goal  of  portability  by  the  foot  soldier  has  been  a 
major  factor. 

The  ASL  AN/GVS-5  based  visioceilometer  incorporates 
the  following  system  design  parameters : 


laser  wavelength: 

1.06 

ym 

laser  pulse  length: 

6  ns 

laser  pulse  energy: 

10  mJ 

source  diameter: 

^1.6 

cm 

source  divergence: 

1.0 

mr 

receiver  diameter : 

5.72 

cm 

receiver  f ield-of-view: 

3.0 

mr 

The  ASL  Model-1  represents  one  alternative  approach 
and  motivation  in  the  development  of  a  prototype  shipboard 
LIDAR-based  visibility  measurement  system.  The  emphasis 
in  the  development  of  this  device  has  been  on  compactness 
and  portability  for  use  by  the  foot-soldier.  The  incor¬ 
poration  of  MIL-SPEC  ruggedized  hardware  components  such 
as  the  ANGVS/5  laser  range  finder  is  attractive  for  field 
or  shipboard  use  since  these  components  are  the  result  of 
extensive  engineering  investments  and  are  designed  to  oper¬ 
ate  under  a  wide  range  of  environmental  conditions  and 
treatment  for  which  laboratory  or  research  grade  equipment 
is  not  suited. 


83 


The  ASL  Model-1  visioceilometer  was  recently  used  in 
conjunction  with  tests  carried  out  at  the  NOSC  remote 
sensing  facility  on  Point  Loma,  San  Diego,  CA  [461  .  The 
limited  comparisons  of  vertical  aerosol  extinction  profiles 
measured  with  an  aircraft-mounted  aerosol  spectrometer  and 
vertical  soundings  with  the  Model-1  LIDAR  show  good  agree¬ 
ment  in  one  case  and  disagreement  in  another,  attributable 
to  possible  differences  in  the  aerosols  along  the  paths 
sampled  by  the  vertically  spiraling  aircraft  and  by  the 
LIDAR.  The  position  taken  by  the  authors  is  that  the 
cursory  comparisons  attempted  thus  far  are  inadequate  to  draw 
meaningful  conclusions  and  that  further  tests  will  be  needed 
to  discern  any  systematic  differences  in  comparisons  of  the 
two  approaches  to  extinction  measurements. 

3.4  DATA  PROCESSING  OF  PULSED  LIDAR  SIGNALS 

3.4.1  CHARGE  COUPLED  DEVICE  (CCD)  -  MICROPROCESSOR 
APPROACH  IN  THE  ASL  MODEL-L  VISIOCEILOMETER 

The  signal  processing  used  by  the  Model-1  visioceilo¬ 
meter  involves  several  distinct  parts.  These  serve  to 
record  the  range  dependent  LIDAR  return  for  either  display 
or  analysis.  The  two  basic  problems  with  recording  a 
LIDAR  return  are  the  large  dynamic  range  of  the  signal  and 
its  wide  bandwidth.  The  first  of  these  problems  is  ad¬ 
dressed  by  a  logarithmic  gain  compression  scheme  which 
compresses  four  decades  of  input  to  two  decades  of  output. 
This  compressed  signal  is  then  recorded  in  a  rather  spe¬ 
cial  way.  Because  the  visioceilometer  is  intended  to  be 
a  hand-held,  battery  powered  unit,  its  power  consumption 
must  be  low.  On  the  other  hand,  the  bandwidth  of  the  sig¬ 
nal  to  be  recorded  must  be  very  large  since  digital  re¬ 
cording  requires  high  sample  rates  (such  as  20  MHz) .  Un¬ 
fortunately,  the  analog  to  digital  (A/D)  conversion  devices 
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which  operate  at  the  required  speeds  and  accuracy  consume 
rather  large  amounts  of  power  and  are  not  well  suited  for 
a  battery  powered,  hand-held  device.  This  apparently  led 
to  the  use  of  a  charge  coupled  device  (CCD)  as  a  high  speed 
analog  buffer.  Analog  samples  of  the  compressed  detector 
video  are  placed  into  the  CCD  buffer  at  a  20  MHz  rate.  The 
CCD  accomplishes  this  with  very  low  power  consumption.  Af¬ 
ter  the  signal  is  recorded  in  this  manner,  it  can  be  read 
out  of  the  CCD  at  a  slower  rate  (44  KHz  in  the  visioceilo- 
meter  case)  such  that  a  low  power  A/D  converter  may  be  used. 
Once  the  digital  version  of  the  recorded  signal  is  obtained, 

it  is  post-processed  by  a  minicomputer  for  display  or  analy- 

2 

sis.  Such  things  as  the  1/r  dependence  of  the  backscat- 
tered  signal  are  removed  at  this  point. 

While  use  of  the  CCD  is  critical  to  the  hand-held  unit, 
it  does  add  some  complications.  A  major  one  is  a  time  de¬ 
pendent  decay  of  the  signal  held  in  the  CCD  buffer.  Each 
element  of  the  signal  is  held  in  the  buffer  for  a  different 
length  of  time.  Thus,  this  decay  appears  as  a  range  depen¬ 
dent  variation  which  must  be  removed  in  post  processing. 
Unfortunately,  this  decay  is  temperature  dependent  such 
that  its  accurate  removal  is  not  trivial. 

For  application  to  a  shipboard  system,  however,  we 
are  not  faced  with  the  same  stringent  power  requirements. 
Thus,  direct  digitization  of  the  compressed  detector  video 
is  a  possibility.  This  would  eliminate  the  temperature 
dependent  decay  problems  associated  with  the  CCD.  Addi¬ 
tionally,  other  signal  processing  techniques  could  be  use- 

2 

ful.  As  an  example,  the  1/r  characteristic  of  the  signal 
could  be  removed  in  real  time  prior  to  the  logarithmic  com¬ 
pression  and  digitization.  This  would  greatly  simplify 
the  required  post  processing  of  the  LIDAR  return. 
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3.4.2  HIGH  SPEED  ANALOG-TO-DIGITAL  CONVERSION  (ADC) 

-  MINICOMPUTER  APPROACH 

The  design  philosophy  which  has  been  emphasized  in  the 
ASL  visioceilometer  data  acquisition  subsystem  design  de¬ 
scribed  in  the  previous  section  has  emphasized  low  power 
consumption  and  miniaturization.  The  use  of  a  CCD  for 
transient  analog  sampling  of  the  LIDAR  return,  followed  by 
slower  analog-to-digital  conversion  (ADC)  is  well  suited  to 
the  configuration  of  that  system  apart  from  the  thermal  sta¬ 
bility  problems  mentioned.  In  a  shipboard  LIDAR  system 
design  the  constraints  of  space  and  power  consumption, 
while  important,  are  nowhere  as  severe  as  in  a  unit  designed 
for  hand-held  operation  as  in  the  visioceilometer. 

High-performance,  compact  ADC  modules  are  currently 
available  which  exceed  the  20  MHz  sampling  rate  used  in 
the  visioceilometer. 

A  modularized  multi-channel  analyzer  system  offered 
by  Le  Croy  Research  Systems  Corporation  (model  3500)  in¬ 
cludes  a  programmable  8085  A  microprocessor,  display,  data 
storage,  and  CAMAC-based  modular  data  acquisition  features. 
Small,  plug-in  transient  digitizer  modules  compatible  with 
this  system  are  currently  available  with  true  8-bit  dynamic 
range  at  50  MHz  sampling  rates.  This  is  over  twice  the 
sampling  rate  currently  used  in  the  ASL  Model-1  visioceilo- 
meter.  A  very  high-speed  ADC  module  for  this  system  will 
be  available  by  summer  1982  offering  a  200  MHz  sampling  rate 
with  true  8-bit  analog  to  digital  conversion  [47] . 

The  large  variations  seen  in  the  data  presented  in 
Section  2.2  and  the  lack  of  any  clear  choice  of  model  to 
represent  these  variations  argue  for  continued  collection 
of  LIDAR  measurements  in  maritime  conditions  and  evaluation 
of  these  data.  To  the  extent  that  the  LIDAR  returns  are 
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processed  on-line  with  pre-selected  algorithms  using  pre¬ 
programmed  microprocessors,  the  value  of  the  data  for  es¬ 
tablishing  the  suitability  of  various  signal  inversion 
methods  (e.g.  the  Klett  inversion  method)  is  reduced.  Cur¬ 
rently  the  ASL  Model-1  visioceilometer  data  are  reduced 
off-line  using  a  minicomputer  system  [46]  but  eventually 
on-line  microprocessor  reduction  of  the  LIDAR  returns  is 
envisioned.  The  latter  approach  will  be  required  in  a 
shipboard  system  designed  for  routine  aerosol  extinction 
measurement.  Based  on  the  current  state  of  understanding 
of  experimental  data  and  model  development,  it  seems  appro 
priate  to  exercise  a  shipboard  LIDAR  system  at  one,  or  pre 
ferably  two  infrared  wavelengths,  and  to  test  the  validity 
of  various  signal  inversion  approaches. 
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4 

LIDAR  APPLICATION  TO  A  SHIPBOARD  MEASUREMENT  ENVIRONMENT 

4 . 1  EYE  SAFETY 

Eye  safety  concerns  and  potential  problems  are  an  im¬ 
portant  aspect  of  the  design  and  operation  of  a  pulsed 
laser  system  especially  when  it  is  intended  for  use  in  an 
outdoor  environment.  A  shipboard  measurement  system  based 
on  the  LIDAR  approach  would  necessarily  have  to  be  operated 
in  close  proximity  to  ship ' s  crew  members  routinely  per¬ 
forming  a  wide  variety  of  tasks  taking  place  at  many  loca¬ 
tions  throughout  the  ship.  To  the  extent  that  it  is  pos¬ 
sible  to  design  an  eye-safe  LIDAR  system  and  retain  adequate 
system  measurement  sensitivity,  a  significant  advantage 
can  be  realized.  Approval  for  fleet  utilization  can  be 
greatly  facilitated  or  impeded  by  the  eye  safety  issue. 

The  apparent  reduction  in  sensitivity  arising  from  reduced 
system  aerosol  backscatter  efficiency  for  wavelengths  longer 
than  1.4  ym  can  be  more  than  offset  by  the  increase  in  maxi¬ 
mum  permissible  exposure  (MPE)  (for  direct  intrabeam  view¬ 
ing)  of  ^  2,000  with  respect  to  the  Nd-YAG  laser  wavelength 
of  1.06  ym. 

Instrument  optical  alignment  procedures  and  adjustments 
which  involve  working  at  close  quarters  with  a  1.06  ym 
laser  source  and  consequently  include  the  possibility  of 
direct  specular  reflection  of  a  small  diameter  beam  into 
the  eye  are  especially  potentially  hazardous.  For  visible 
sources  with  wavelengths  below  0.7  ym  the  operator's  eye 
can  sense  and  help  in  avoidance  of  direct  viewing  of  the 
beam.  However  coherent  laser  radiations  at  wavelengths 
between  0.7  ym  and  1.4  ym  do  not  cause  a  visual  sensation 
in  spite  of  being  directly  transmitted  by  the  lens  of  the 
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eye  and  imaged  as  a  potentially  damaging  intense  spot  on 
the  retina,  requiring  that  particular  caution  must  be  exer¬ 
cised  while  working  with  such  sources.  Evaluation  of  al¬ 
ternatives  to  the  use  of  readily  available  and  reliable  Nd- 
YAG  laser  sources  should  be  given  careful  consideration 
from  the  eye  safety  viewpoint  alone,  apart  from  the  many 
other  issues  to  be  considered  in  the  development  of  a  ship¬ 
board  measurement  system. 

4.1.1  ESTIMATE  OF  EYE  EXPOSURE  HAZARD  FROM  THE  ASL 
MODEL-1  VISIOCEILOMETER 

The  maximum  permissible  exposure  (MPE)  for  intrabeam 

viewing  of  a  1.06  ym  pulsed  laser  (single  pulse  between 

•9  -5  2 

10  and  5  x  10  sec)  is  5  x  10  J/cm  according  to 

guidance  published  by  the  American  National  Standards  In¬ 
stitute  [48] .  If  the  pulse  repetition  frequency  (PRF)  is 

increased  to  10  pps  then  the  MPE  value  decreases  to 
_6  2 

1.7  x  10  J/cm  per  pulse.  Assuming  a  PRF  of  1  pps  for 
the  AN/GVS-5  specifications  given  in  Section  3.3  then  a 
safe  eye  exposure  distance  (SEED)  may  be  calculated.  The 
beam  divergence  is  specified  as  1.0  mr  which  is  about  6 
times  the  diffraction  limit  of  a  Gaussian  mode  diffracting 
from  the  0.8  cm  aperture .  The  SEED  value  S  can  be  found 
from  the  relation: 


w(f) 


MPE 


(17) 


where  E^  is  the  laser  pulse  energy  in  joules  and  8  i 3  the 
beam  divergence  angle  in  radians.  Using  the  ASL  Model-1 
values  we  find: 


4  x  10 
- - 

ir(S  x  10  -V 


<  5  x  10 


J/cm 


(18) 
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or  S  >  .505  x  10^  cm  =  505  m.  If  the  PRF  is  increased  to 
10  pps  the  SEED  value  increases  to  about  1.6  km.  The  ANSI 
requirement  [48]  states  that  direct  intrabeam  viewing  of 
the  launched  laser  pulse  must  be  prevented  until  the  beam 
has  traveled  at  least  505  m  for  the  1  pps  case  and  1.6  km 
for  the  10  pps  case.  An  alternate  approach  to  this  poten¬ 
tial  eye  safety  problem  can  be  evaluated  by  considering  what 

increase  in  transmitter  aperture  size  would  be  needed  to 

_2 

reduce  the  SEED  value  to  zero  for  a  pulse  energy  of  10  J. 
That  is,  how  large  must  the  launched  beam  cross  sectional 
area  (Ag)  be  before  it  is  eye-safe  as  it  leaves  the  trans¬ 
mitting  aperture?  In  this  case: 


10  3°ules  <  5  x  1(T6  J/cm2  (19) 

a 

3  2 

or  Ag  >  2  x  10  cm  .  This  requires  a  transmitter  mirror 
diameter  of  about  50  cm  (20")  diameter  which  is  clearly  out 
of  the  question  for  the  Army  application  but  could  be  con¬ 
sidered  for  a  shipboard  measurement,  keeping  in  mind  the 
close  quarters  that  such  a  unit  would  be  used  in  aboard 
ship. 


4.1.2  WAVELENGTH  CONTROL  FOR  EYE  SAFETY 

If  an  operating  wavelength  longer  than  1.4  ym  can  be 

used,  the  MPE  value  increases  dramatically  from  5  x  10-® 
2-2  2 

J/cm  to  10  J/cm  for  1  pps  repetition  rates.  For  wave¬ 
lengths  of  1.4  ym  and  longer  the  lens  of  the  human  eye  is 
no  longer  transparent  and  focusing  of  an  exceedingly  small 
and  intense  spot  of  coherent  light  on  the  retina  of  the  eye 
no  longer  occurs.  Possible  eye  hazard  considerations  then 
concern  exposure  of  the  sensitive  surface  of  the  eye  to  the 
laser  radiation  and  depend  upon  the  actual  energy  density 
in  the  beam  at  the  surface  of  the  eye,  not  in  a  tiny  focal 
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spot  formed  on  the  retina.  Consequently  an  increase  of  about 
2000  in  the  MPE  values  applies. 

It  can  be  seen  from  Equation  19  that  a  transmitter 

2 

aperture  area  only  slightly  larger  than  1  cm  would  corres¬ 
pond  to  an  intra-beam  energy  density  meeting  the  MPE  re¬ 
quirement  for  wavelengths  longer  than  1.4  urn  at  1  pps  and 

2 

that  a  transmitted  beam  area  of  about  1.6  cm  would  be  eye 
safe  for  10  pps  repetition  rates. 

4. 1.2.1  Stimulated  Raman  Scattering  as  a  Wavelength 
Shifting  Mechanism 

We  have  seen  from  the  discussion  and  MPE  values  pre¬ 
sented  in  the  previous  section  that  a  significant  reduction 
in  potential  eye  safety  hazards  results  if  laser  wavelengths 
longer  than  1.4  vim  can  be  used  with  a  LIDAR  system.  Several 
examples  of  highly  developed,  MIL-SPEC  laser  hardware  oper¬ 
ating  at  the  Nd-YAG  wavelength  of  1.06  ym  are  currently 
being  widely  used  in  military  laser  designators  and  range¬ 
finders.  These  devices  have  resulted  from  extensive  develop¬ 
ment  programs  designed  to  produce  extremely  rugged  and  re¬ 
liable  components  which  will  operate  satisfactorily  under 
very  adverse  circumstances  and  extremes  of  environmental 
conditions.  The  ASL  Model-1  visioceilometer  is  based  on 
the  use  of  the  AN/GVS-5  laser  rangefinder  source  and  re¬ 
ceiver  modules  which  are  MIL-SPEC  components.  Additional 
MIL-SPEC  systems  will  be  discussed  in  the  following  section. 

It  would  be  highly  advantageous  in  the  development  of 
a  shipboard  LIDAR  to  capitalize  on  the  extensive  hardware 
development  which  has  occurred  in  the  design  and  production 
of  reliable  MIL-SPEC  laser  sources.  Nearly  all  of  these 
systems  operate  at  1.06  ym,  so  that  an  efficient  and  reli¬ 
able  means  of  shifting  the  1.06  ym  wavelength  to  1.4  ym  or 
longer  is  needed  in  order  to  utilize  exsisting  MIL-SPEC 
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hardware  and  to  simultaneously  obtain  an  "eye  safe"  output 
for  use  as  a  LIDAR  source. 

Recent  developments  in  laser  frequency  conversion  us¬ 
ing  Stimulated  Raman  Scattering  (SRS)  have  shown  that  energy 
conversion  efficiencies  to  the  first-Stokes  output  wave¬ 
lengths  for  high  pressure  gaseous  D2  and  CH4  exceed  40%  [49] . 
Design  considerations  and  test  results  for  a  multiple-pass, 
optical-resonator  Raman  cell  have  been  presented  recently 
by  Trutna  and  Byer  [50] .  Application  of  these  techniques 
to  the  design  and  testing  of  a  compact,  efficient  SRS  cell 
used  to  shift  1.06  ym  laser  radiation  to  1.54  ym  using  a 
methane  filled  cell  and  to  1.9  ym  using  hydrogen  has  been 
carried  out  by  Hughes  Research  Laboratories  [51] .  The  SRS 
first-Stokes- shifted  wavelength  for  CH.  pumped  with  the 
output  of  a  Nd-YAG  laser  is  1.543  ym  (6480.68  cm  ).  If  a 
deuterium  cell  were  used  the  output  wavelength  would  be 
approximately  1.561  ym  (6407.88  cm-1).  The  SRS  wavelength 
shifts  are  fixed  by  the  vibrational  energies  of  the  scat¬ 
tering  molecules  (2914.2  cm-1  for  CH4  and  2987  cm-1  for  D2) . 
Slight  tuning  adjustments  of  the  SRS  output  frequency  can 
be  made  to  avoid  possible  overlap  with  interfering  atmos¬ 
pheric  absorption  lines  by  use  of  an  etalon  and/or  tempera¬ 
ture  tuning  the  1.064  ym  Nd-YAG  output  wavelength  over  a 
few  tenths  of  a  wavenumber. 

Figures  30  and  31  show  the  calculated  atmospheric 
transmission  of  a  5  km  sea  level  path  for  a  mid-latitude 
summer  atmosphere  (14.25  torr  partial  pressure  of  H20)  in 
the  spectral  region  near  1.54  ym.  The  interval  between 
6460  cm  1  and  6500  cm  1  is  shown  in  Figure  30;  Figure  31 
is  a  higher  dispersion  plot  of  the  interval  between  6473  cm"1 
and  6483  cm  1  showing  the  location  of  the  1.543  ym  first- 
Stokes  SRS  line  for  CH^ .  From  an  examination  of  Figure  31 
it  can  be  seen  that  a  prominent  atmospheric  absorption  feature 
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occurs  near  the  1.543  ym  SRS  line.  The  major  contribution 

to  this  absorption  feature  is  an  H-0  line  located  at 

6481.1110  cm  ;  minor  contributions  are  due  to  weaker  CO., 

-1  -l 

and  ^0  lines  located  at  6481.366  cm  and  6481.564  cm 

respectively  [52] .  By  tuning  the  Nd-YAG/CH.  SRS  output  line 

-1  -1  .  4 
by  0.2  cm  to  0 . 4  cm  ,  this  near-coincidence  of  the  CH4 

SRS  line  and  the  atmospheric  water  vapor  line  could  be  used 
as  the  basis  of  a  differential  absorption  LIDAR  (DIAL)  scheme, 
expanding  the  usefulness  of  the  shipboard  LIDAR  to  include 
a  capability  for  remote,  long-path  measurements  of  water 
vapor  density  in  addition  to  the  primary  measurement  function 
of  aerosol  extinction.  It  should  be  feasible  (e.g.  by  means 
of  an  oscillating  etalon)  to  shift  the  SRS  line  on  and  off 
of  the  atmospheric  H20  line  on  alternate  pulses  thereby 
providing  two  distinct  LIDAR  returns  which  can  be  ratioed 
to  give  the  water  vapor  attenuation  profile.  Whether  or 
not  this  particular  aspect  of  the  1.543  ym  SRS  line  of  the 
Nd-YAG/CH4  system  is  exploited,  the  SRS  technique  shows  a 
great  deal  of  promise  in  providing  a  relatively  straight¬ 
forward  means  to  generate  "eye  safe"  laser  radiation  using 
a  relatively  passive  Raman  cell  together  with  proven,  high¬ 
ly  developed,  MIL-SPEC  1.06  ym  laser  hardware. 

4.2  SUITABILITY  OF  EXSISTING  MIL-SPEC  AND  OTHER  RUGGEDIZED 

E-0  HARDWARE 

4.2.1  MIL-SPEC  1.06  ym  SYSTEMS 

The  AN/GVS-5  laser  rangefinder  system  which  has  been 
incorporated  into  the  Army  ASL  Model-1  visioceilometer  is 
an  example  of  the  use  of  MIL-SPEC  components  for  a  LIDAR 
system.  The  characteristics  of  the  AN/GVS-5  were  listed 
in  Section  3.3.  Other,  higher  energy  1.06  ym  MIL-SPEC 
sources  have  been  developed  for  military  laser  designator 
applications  and  could  be  utilized  in  the  development  of  a 
shipboard  LIDAR.  One  such  device  is  the  AN/TVQ-2  Army 
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Ground  Laser  Locator  Designator  (GLLD) .  The  GLLD  output 
energy  is  about  a  factor  of  ten  larger  than  that  of  the 
AN/GVS-5  and  the  device  characteristics  include  reasonably 
good  beam  quality.  The  combined  weight  of  the  laser  desig¬ 
nator-rangefinder  module  of  the  GLLD  is  12.7  kg  (28  lbs.) 
and  the  battery  supply  module  providing  completely  self- 
contained  operation  weighs  3.2  kg  (7  lbs.).  While  this 
device  is  substantially  larger  than  the  AN/GVS-5,  its'  size 
and  weight  are  quite  compatible  with  use  in  a  shipboard 
operation.  The  most  probabl  .  configuration  for  using  it 
as  part  of  a  shipboard  LIDAR  would  involve  interfacing  the 
GLLD  laser  source  with  moderately  large  (^30  cm)  aperture 
transmit/receive  telescope  optics.  A  more  detailed  discus¬ 
sion  of  such  a  configuration  together  with  performance  es¬ 
timates  are  presented  in  Section  4.3.  Use  of  an  SRS  cell 
together  with  the  GLLD  source  to  provide  output  pulse  ener- 
geries  in  the  range  of  tens  of  milli joules  at  1.54  urn  offers 
several  attractive  advantages.  Coolina  for  the  GLLD  laser 
is  provided  by  pressurized  N2  gas  circulated  through  a  N2~ 
air  heat  exchanger  through  which  filtered  external  air  is 
circulated  by  means  of  a  blower.  This  technique  offers  the 
advantage  that  no  cooling  liquids  are  needed.  The  entire 
laser  optical  train  and  N2  circulation  path  are  hermetically 
sealed  against  the  intrusion  of  contamination  which  could 
rapidly  render  the  laser  system  inoperative.  Thus,  these 
features  of  the  device  are  ideally  suited  to  a  shipboard 
installation  where  high  humidity  and  salt-laden  condensates 
can  prove  very  detrimental  to  the  delicate  and  critical 
optical  components  comprising  the  laser  optical  cavity. 

Another  MIL-SPEC  1.06  urn  laser  designator  system  in¬ 
corporating  similar  design  features  to  the  GLLD  but  with 
about  twice  as  large  an  output  pulse  energy  is  used  in  the 
Navy  A- 6  aircraft  TRAM  system.  The  N2~air  heat  exchanger 
feature  is  used  in  the  TRAM  designator  as  well. 
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Both  the  GLLD  system  and  the  TRAM  system  are  fully 
developed  and  currently  in  or  near  production.  Availability 
of  engineering  development  units  from  either  program  for 
evaluation  as  to  their  potential  for  use  in  a  shipboard 
LIDAR  development  program  should  be  pursued. 

Several  other  liquid-cooled  MIL-SPEC  1.06  urn  sources 
are  in  production  or  advanced  development  and  warrant  fur¬ 
ther  investigation.  At  this  time  the  absence  of  a  liquid 
coolant  and  the  compactness  and  ruggedness  of  the  GLLD  and 
TRAM  designator  system  designs  appear  to  offer  the  most  de- 
sireable  characteristics  for  use  in  a  shipboard  LIDAR. 

4.2.2  SYSTEMS  OPERATING  AT  10.6  pm 

An  alternative  to  operation  at  1.06  ym  or  at  the  SRS 
shifted  wavelength  of  1.54  ym  is  the  use  of  a  CC>2  laser 
operating  at  10.6  ym.  This  wavelength  is  well  into  the 
eye  safe  region  but  the  anticipated  LIDAR  backscatter  re¬ 
turn  signal  is  expected  to  be  substantially  smaller  in  most 
cases  than  the  return  signal  at  1.06  ym  or  1.54  ym  (see 
Section  2.4,  Figures  27  and  28).  Cryogenically  cooled  de¬ 
tectors  will  be  required  with  a  10.6  ym  system  and  hetero¬ 
dyne  detector  techniques  can  be  employed  to  gain  a  few 
orders  of  magnitude  in  system  sensitivity  to  offset  the 
low  levels  of  aerosol  backscatter  at  10.6  ym. 

Laser  pulse  energies  greater  than  10  1  joule  are  de- 
sireable  for  operation  in  moderate  visibility  conditions 
out  to  a  few  kilometers.  At  the  present  time  development 
of  MIL-SPEC  C02  laser  sources  has  lagged  substantially  be¬ 
hind  that  of  1.06  ym  sources  for  US  military  systems.  A 
promising,  albeit  rather  low  energy,  ruggedized  C02,  trans¬ 
verse-electrical-atmospheric  pressure  (TEA)  system  is  manu¬ 
factured  in  the  UK.  This  laser  is  used  in  the  Marconi 
Avionics  and  Ferranti  model  307  C02  laser  rangefinders 
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manufactured  for  the  military.  The  laser  is  of  a  rugged- 
ized  design,  well  suited  for  field  use.  The  output  pulse 
energy,  however  is  specified  at  about  13  x  10~3J  with  a 
peak  power  of  220  kw  and  pulse  duration  of  60  ns,  making  its 
use  as  a  source  for  10.6  urn  aerosol  backscatter  measurements 
marginal  for  all  but  low  visibility  fog-type  conditions 
where  scattering  at  10.6  ym  is  more  nearly  comparable  to 
that  at  1.06  ym. 

Several  developmental ,  compact  C02  TEA  lasers  are  being 
used  in  LIDAR  type  research  instruments.  These  sources 
typically  operate  at  higher  pulse  energy  than  the  UK  source 
described  above  but  are  not  as  well  suited  for  routine  use 
in  the  field.  The  most  promising  approach  at  the  present 
time  for  obtaining  a  10.6  ym  source  suitable  for  use  in 
a  shipboard  LIDAR  would  seem  to  involve  the  upgrade  and 
environmental  protection  of  a  compact,  energetic  TEA  laser 
system  operating  in  10  1  to  1.0  J/pulse  energy  range. 

4.3  SHIPBOARD  LIDAR  SYSTEM  PERFORMANCE  ESTIMATES 

4.3.1  ASL  MODEL-1  VISIOCEILOMETER  SYSTEM 

Using  the  results  of  calculations  discussed  earlier 
in  Section  2.4  and  shown  in  Figure  32  together  with  the 
ASL  Model-1  visioceilometer  characteristics  discussed  in 
Section  3.3,  performance  estimates  of  this  system  were 
carried  out.  The  results  of  these  estimates  are  shown  in 
Figure  32  for  ranges  between  200  m  and  1.0  km.  Calculated 
peak  power  of  the  LIDAR  return  signal  in  watts  is  shown 
as  a  function  of  range  for  four  model  maritime  aerosol  dis¬ 
tributions  with  2  km,  5  km  and  10  km  visibilities  and  70% 
relative  humidity.  Also  shown  by  the  dashed  curve  is  the 
energy  density  in  the  transmitted  beam  E '  as  a  function  of 

range.  One  can  see  that  E'  is  substantially  larger  than 

—6  2 

the  1.06  ym  MPE  value  of  5  x  10  J/cm  (single  pulse)  or 
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32.  LIDAR  TRANSMITTED  ENERGY  DENSITY  E' 
AND  RECEIVED  PEAK  POWER  PD  VERSUS 
RANGE  FOR  THE  ASL  MODEL-1  (AN/GVS-5) 
VISIOCEILOMETER  SYSTEM  FOR  VARIOUS 
TYPES  OF  MARINE  AEROSOLS. 


(j/CM 


“6  2 

1.6  x  10  J/cm  at  10  pps  for  ranges  shorter  than  about 

500m.  The  safe  eye  exposure  distance  (SEED)  for  this  device 

was  calculated  in  Section  4.1.1  and  found  to  be  505  m  for 

1  pps  and  1.6  km  for  10  pps.  The  return  signal  peak  powers 

are  seen  to  vary  between  10  7  watts  at  'v  2  00  m  to  less  than 
-9 

10  watts  for  ranges  in  excess  of  1  km,  indicating  that 
system  performance  will  be  marginal  for  greater  ranges  in 
intermediate  backscatter  situations. 

4.3.2  UPGRADED  1.06  urn  SYSTEM 

The  same  type  of  analysis  as  shown  in  "’igure  32  for 
the  ASL  Model-1  system  was  carried  out  for  a  hypothetical 
upgraded  system  with  the  following  characteristics:  pulse 
energy,  Eq  =  10-^  J/pulse  and  transmitter  and  receiver 
diameter  of  30.5  cm  and  a  beam  divergence  of  1.0  mrad,  the 
same  as  the  ASL  Model-1  system.  The  results  of  these  cal¬ 
culations  are  shown  in  Figure  33.  The  larger  transmitter 
aperture  more  than  offsets  the  factor  of  ten  pulse  energy 
increase  to  give  a  lower  energy  density  at  the  device  out¬ 
put  aperture,  however  it  is  still  not  below  the  eye  safe 
—  6  2 

limit  of  5  x  10  J/cm  .  As  can  be  seen  in  Figure  33  the 
energy  density  in  the  beam  for  this  system  is  still  deter¬ 
mined  primarily  by  the  beam  divergence  and  remains  above 
the  eye  safe  level  over  the  ranges  (0.2  km  to  1.0  km)  shown 

in  the  figure.  If  the  sytem  operating  wavelength  were 

-2  2 

longer  than  1.4  ym  the  MPE  would  be  10  J/cm  and  there 
would  be  no  eye  hazard  even  at  the  device  aperture.  The 
LIDAR  return  signal  levels  for  the  upgraded  configuration 
shown  in  Figure  33  are  approximately  2  orders  of  magnitude 
greater  than  those  for  the  baseline  ASL  Model-1  device  in¬ 
dicating  that  substantially  improved  performance  can  be 
expected  for  ranges  longer  than  1  km  when  compared  to  the 
ASL  Model-1  system. 
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FIGURE  33.  LIDAR  TRANSMITTED  ENERGY  DENSITY  E' 
AND  RECEIVED  PEAK  POWER  PD  VERSUS 
RANGE  FOR  AN  UPGRADED  LIDAR  SYSTEM 
WITH  THE  CHARACTERISTICS  SHOWN  IN 
THE  FIGURE. 
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Figure  34  shows  a  performance  comparison  of  the  two 
systems  for  fog  and  2  km  visibility  oceanic  aerosol  condi¬ 
tions  for  ranges  between  500  m  and  5  km.  The  upgraded 
system  is  seen  to  generate  useful  return  signal  levels  out 
to  5  km  for  the  oceanic  aerosol  while  the  AN/GVS-5  based 
system  (ASL  Model-1  visioceilometer )  is  seen  to  generate 
marginally  useful  return  signal  levels  for  ranges  greater 
than  about  1.5  km . 

4.3.3  COMPARISON  OF  SYSTEMS  OPERATING  IN  THE  1.06  ym, 
3-5  ym,  and  8-12  ym  REGIONS 

Using  the  upgraded  system  characteristics  and  the  cal¬ 
culated  backscatter  coefficients  for  a  2  km  visibility 
oceanic  aerosol  a  system  performance  comparison  for  3 
different  wavelengths  was  generated.  The  results  of  this 
comparison  are  shown  in  Figure  35  for  ranges  between  200  m 
and  1.0  km.  The  rapid  decrease  of  return  signal  peak  power 
with  increasing  wavelength  shows  that  higher  output  pulse 
energies  and  increased  detector  sensitivity  e.g.  a  hetero¬ 
dyne  configuration  will  be  required  for  systems  operating 
at  10.6  ym. 

4.3.4  DEPENDENCE  OF  LIDAR  SYSTEM  PERFORMANCE  ESTI¬ 
MATES  UPON  AEROSOL  EXTINCTION  DATA  AND  MODEL 
VARIABILITY 

When  the  variability  in  the  infrared  aerosol  extinc¬ 
tion  data  described  in  Section  2.2  is  taken  into  account , 
adjustments  to  the  nominal  case  backscatter  calculations 
shown  in  Figure  27  can  be  made.  Based  on  the  range  of  values 
for  a  seen  in  the  experimental  aerosol  extinction  data  one 
could  anticipate  worst  case  corrections  to  the  calculated 
data  shown  in  Figure  27  of  about  0.1  at  a  wavelength  of 
0.55  ym,  0.08  at  1.06  ym  and  0.25  at  3.8  ym.  Lacking  a 
sufficient  body  of  10.6  ym  maritime  optical  transmission 
data  and  considering  the  uncertainties  associated  with  the 
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FIGURE  34.  COMPARISON  OF  RECEIVED  PEAK  POWER  PD 
FOR  THE  AN/GVS-5  BASED  AND  UPGRADED 
SYSTEMS  VERSUS  RANGE. 


RANGE  (KM) 


FIGURE  35.  COMPARISON  OF  RECEIVED  PEAK  POWER  LEVELS; 

PD  VERSUS  RANGE  FOR  THREE  SYSTEMS  OPERATING 
AT  PULSE  ENERGIES  OF  10"1  J/PULSE  WITH  30.5 
cm  DIAMETER  COLLECTING  OPTICS  AND  1.06  ym, 
3.8  ym  AND  10.6  ym  WAVELENGTHS . 


modeling  of  the  water  vapor  continuum  absorption  in  the 
8-12  ym  spectral  region  [23],  the  0.25  correction  factor 
appropriate  for  3.8  ym  may  be  applied  to  system  performance 
estimates  for  10.6  ym  as  well. 

Model  predictions  for  ajR/aVIS  have  been  seen  to  vary 
by  about  a  factor  of  two  with  changes  in  relative  humidity 
and  by  about  a  factor  of  four  with  changes  in  windspeed. 
Considering  the  observed  variations  in  experimental  data 
and  the  variations  exhibited  by  the  computed  model  values 
for  0IR  over  the  range  of  windspeeds  and  relative  humidi¬ 
ties  of  greatest  interest,  the  performance  predictions 
shown  in  Figures  32-35  should  reliable  at  the  specified 
levels  or  somewhat  lower  but  by  no  more  than  a  factor  of 
ten. 


4.3.5  SIGNAL  TO  NOISE  RATIO  ESTIMATES 

Using  the  estimates  of  collected  signal  power  for 
various  system  configurations  and  ranges  presented  in 
Figures  32-35  together  with  detector  noise  equivalent 
power (NEP)  estimates,  the  signal  to  noise  ratio  achievable 
with  a  prototype  LIDAR  system  can  be  estimated.  Table  7 
summarizes  signal  to  noise  ratio  calculations  for  several 
of  the  shipboard  LIDAR  configurations  previously  discussed. 
Column  one  of  the  table  lists  the  system  configuration. 

The  ASL-1  system  is  the  visioceilometer  device  based  on  the 
AN/GVS-5  laser  rangefinder  which  is  described  in  section  3.3. 
The  upgraded  system  configuration  listed  in  column  one  is 
based  on  a  transmitter  pulse  energy  of  10-1  joule  and  a 
receiver  collector  diameter  of  30.5  cm.  Signal  to  noise 
ratio  estimates  for  several  system  configurations  including 
upgraded  1.06  ym,  1.54  ym,  3.8  ym  and  10.6  ym  direct  and 
heterodyne  detection  systems  are  shown  in  the  table.  Column 
two  lists  the  operating  wavelength.  Column  three  lists  the 
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TABLE  7 

SIGNAL  TO  NOISE  RATIO 
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appropriate  detector  and  its  operating  temperature  is 
listed  in  column  four.  The  detector  D*,  operating  bandwidth 
and  noise  equivalent  power  are  listed  in  the  next  three 
columns.  Calculated  signal  levels  and  signal  to  noise 
ratios  are  listed  in  the  next  three  sections  of  the  table 
for  three  ranges:  200m,  1  km  and  5  km. 

All  of  the  calculations  are  appropriate  for  a  measure¬ 
ment  bandwidth  of  100  MHz  except  for  the  ASL  model-1  which 
employs  a  50  MHz  bandwidth.  One  can  see  from  an  inspection 
of  the  table  that  the  detector  detectivity  (D*)  is  a  factor 
limiting  performance  at  the  eye  safe  wavelength  of  1.54  pm. 

It  should  be  noted  that  rapid  development  of  detectors  for 
use  at  this  wavelength  is  curren'-ly  taking  place  accellerated 
by  requirements  for  fiber  optical  communication  applications 
(most  optical  fibers  exhibit  very  low  loss  at  this  wavelength) 
and  therefore  that  the  D*  value  contained  in  Table  6  for  the 
Ge  photodiode  detector  will  very  likely  be  exceeded  in  the 
near  future  by  some  of  the  newer  InGaAs  photodiodes  being 
optimized  for  use  at  this  wavelength. 

One  can  see  from  an  inspection  of  the  table  that  all 
of  the  system  configurations  provide  useable  signal  to  noise 
ratios  for  a  1  km  range  with  the  exception  of  the  10.6  pm 
direct  detection  configuration.  For  a  range  of  5  km  most  of 
the  systems  generate  marginally  useful  return  signal  levels 
and  consequently  exhibit  small  signal  to  noise  ratios.  The 
best  performance  is  seen  in  the  upgraded  1.06  pm  system  and 
the  heterodyne  10.6  pm  system.  In  order  to  increase  signal 
to  noise  ratios  co-addition  of  multiple  LIDAR  signal  returns 
can  be  implemented.  The  signal  to  noise  ratio  (S/N)  is  ex- 
pected  to  improve  as  n  '  so  that  S/N  ratios  of  10  should 
be  considered  useful  since  for  n=100,  multiple  laser  pulse 
returns  can  be  easily  coadded  using  a  1  pps  or  10  pps  repetition 
rate.  However  when  10^  or  more  pulses  must  be  coadded  to  achieve 
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an  improved  S/N  of  ^1.0  and  the  maximum  pulse  repetition  rates 
for  the  10  ^  J/pulse  lasers  oeing  considered  are  about  10  pps , 
the  required  co-addition  times  become  prohibitively  long. 
Accordingly,  the  upgraded  1.06  pm  and  3.8  um  configurations 
appear  to  be  useful  for  5  km  ranges.  The  ASL  model  -  1 
system  will  probably  not  be  routinely  useable  at  this  range 
nor  will  the  1.54  pm  upgraded  system  and  the  10.6  um  direct 
detection  system  unless  appreciably  higher  D*  detectors  for 
these  wavelengths  can  be  obtained.  The  advantages  of  hetero¬ 
dyne  detection  become  obvious  when  one  examines  the  last 
row  of  the  table.  Although  the  backscattered  return  signals 
at  10.6  pm  are  2  or  3  orders  of  magnitude  smaller  than  those 
for  the  shorter  wavelengths,  the  gain  in  detector  effective 
noise-equivalent-power  (NEP)  more  than  offsets  this  loss 
and  the  result  is  useable  S/N  ratios  out  to  5  km  for  the 
conditions  indicated  in  Table  6.  The  obvious  conclusions 
that  can  be  drawn  from  the  values  listed  in  the  table  are 
that  improved  detectors  at  1.54  pm  (with  respect  to  the  D* 
value  listed)  and  heteorodyne  detection  at  10.6  pm  will  be 
very  crucial  to  the  successful  implementation  of  LIDAR 
systems  operating  in  the  10  ^  J/pulse  range  at  these  wave¬ 
lengths. 
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4. 3. 5.1  Techniques  for  Signal  to  Noise  Ratio  Improvement 

The  development  of  improved  detectors  for  1.54  um  should 
be  accellerated  by  the  rapid  growth  of  the  fiber  optics 
industry.  Improved  Ge  FIN  photodiodes  and/or  InGaAs 
photodiodes  with  quantum  efficiencies  around  60%  are  being 
developed  and  becoming  available  for  use  in  this  wavelength 
region,  therefore  it  should  be  possible  to  substantially 
improve  S/N  ratios  for  the  1.54  ym  upgraded  system  con¬ 
figuration  shown  in  Table  6.  An  improvement  of  a  factor 
of  ten  in  detector  D*  will  be  needed  for  useful  performance 
of  this  system  at  a  range  of  5  Jon,  increasing  the  estimated 
single  pulse  S/N  to  about  5  x  10  .  LIDAR  return  S/N  ratio 

can  be  further  improved  to  a  useful  level  for  measurement 
by  the  co-addition  of  multiple  pulses.  Assuming  a  10  pps 
repetition  rate,  co-addition  of  return  signals  for  an 
interval  of  40  seconds  (i.e.  400  return  pulses)  would  provide 
a  S/N  ratio  improvement  of  20,  resulting  in  a  final  S/N 
ratio  of  unity. 

There  are  very  large  potential  advantages  to  be  gained 

by  the  use  of  heterodyne  detection  techniques  as  shown  by 

-12 

the  entries  listed  in  the  last  row  of  Table  6  (NEP  ^10  w 

g 

in  a  10  Hz  bandwidth) .  Using  the  heterodyne  approach  the 
backscattered  return  signal  from  the  transmitted  pulse  is 
mixed  with  the  output  of  a  local  oscillator  laser  offset  in 
frequency  by  10  GHz  or  more.  The  LIDAR  return  signal  is 
then  detected  as  the  beat  frequency  difference  signal.  Very 
high  speed  photodetectors  and  amplifiers  are  required  in 
this  configuration.  Use  of  a  rare  isotope  CO2  laser  and 
selection  of  the  respective  laser  transitions  such  that  the 
transmitter  and  local  oscillator  frequencies  are  offset  by 
about  10  GHz  should  be  relatively  straightforward.  Most  of 
the  rare  isotope  COj  laser  lines  have  been  measured  with 
sufficient  accuracy  [53]  so  that  choice  of  appropriate  pairs 
of  operating  lines  should  not  be  a  problem. 
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The  effects  of  laser  speckle  patterns  superimposed  on 
the  spatial  distribution  of  backscattered  radiation  will 
tend  to  deteriorate  its  coherence  and  consequently  will 
partially  ameliorate  the  potential  advantages  offered  by 
heterodyne  detection.  Atmospheric  turbulence  effects  will 
also  cause  some  loss  of  phase  coherence  and  therefore  impact 
heterodyne  receiver  performance.  The  turbulence  effects 
should  tend,  however,  to  contribute  to  some  randomizing  or 
averaging  of  speckle  patterns.  Since  the  two  phenomena 
should  be  uncorrellated  or  possibly  anti-correllated,  the 
resultant  S/N  decrease  due  to  a  noise  like  contribution  by 
laser  speckle,  ngp  and  that  due  to  atmospheric  turbulence, 
n. ,  would  be  expected  to  be  proportional  to  the  root  sum 
square  addition  of  the  two  contributions,  i.e.  n=(n^p  +  n£b)  • 
Some,  if  not  the  full  theoretical  improvement  over  a  direct 
detection  system  should  be  possible. 

-2  2 

The  MPE  for  A>  1.4  ym  is  10  J/cm  for  1  pps.  For  a 
system  with  a  transmitter  aperture  diameter  of  30.5  cm 

2 

the  transmitted  beam  cross  sectional  area  is  about  730  cm  , 
therefore  the  laser  pulse  energy  coul  ;  be  a.»  large  as 
7.3  joules  (for  a  unifo/.mly  .Muz  ;.nated  aperture)  and  the 
system  would  remain  eye  safe.  Clearly  an  increase  in  : 
to  the  1-5  J  range  would  be  permissible.  The  limitations 
in  useable  source  energy  are  seen  to  be  determined  more  by 
available  source  technology  than  by  eye  safety  considerations. 

A  stable  oscillator-TEA  laser  amplifier  configuration  appears 
to  be  the  most  promising  approach  for  obtaining  the  required 
10.6  ym  pulse  energy  at  the  present  time.  Research  LIDAR  systems 
incorporating  these  features  are  currently  being  fielded, 
however  a  significant  amount  of  development  and  packaging 
will  be  needed  before  such  a  system  can  begin  to  approach 
the  reliability  of  the  MIL-SPEC  Nd-YAG  systems  previously 
discussed. 
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4 . 4  SHIPBOARD  IMPLEMENTATION 


Several  of  the  considerations  impacting  the  implemen¬ 
tation  of  a  shipboard  LIDAR  system  have  been  discussed  pre¬ 
viously  in  sections  4.1  and  4.2.  A  summary  of  that  information 
will  be  incorporated  here  by  describing  the  characteristics 
of  a  prototype  shipboard  LIDAR  system,  the  development  of 
which  appears  to  be  both  feasible  and  highly  useful  at  the 
present  time. 

Based  on  the  earlier  discussions  concerning  eye  safety 
and  recognizing  the  ambiguities  present  in  scaling  short 
wavelength  data  to  the  8-12  um  region,  a  system  simultaneously 
operating  at  1.54  um  and  10.6  pm  appears  most  worthwhile. 

Such  a  system  can  be  made  relatively  compact  although  the 
extreme  miniturization  emphasized  in  the  development  of  the 
ASL  model  -  1  system  is  not  required.  A  system  employing 
a  common  30  cm  (12  inch)  diameter  reflecting  telescope  both 
as  transmitter  and  receiver  for  both  wavelengths  appears 
reasonable.  The  optics  and  source  module  for  a  prototype 
shipboard  system  would  be  contained  in  a  cylinder  less  than 
or  equal  to  40  cm  (16  inches)  in  diameter  and  about  1.2  m 
in  length.  A  weight  allowance  of  50  to  75  lbs.  should  be 
sufficient  for  this  assembly.  The  self-contained  data 
acquisition  and  processing  equipment  might  occupy  about 
eight  cubic  feet  and  weigh  about  150  lbs.  This  sub-system 
could  be  remotely  located  from  the  LIDAR  optical  head  con¬ 
taining  optics,  laser  sources,  detectors  and  analog-to- 
digital  converter,  if  desirable.  In  this  configuration 
non-weatherized  electronic  modules  could  be  housed  in  a 
sheltered  area  permitting  system  operation  during  a  wider 
range  of  weather  conditions  than  would  be  possible  if  the 
electronics  sub-system  were  mounted  adjacent  to  the  LIDAR 
optical  head.  Initially  shipboard  testing  of  the  LIDAR 
system  might  best  be  accomplished  using  a  single  physical 
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installation  with  the  data  sub-system  adjacent  to  the 
optical  head  to  permit  simplified  operation  of  the  system. 
Once  an  adequate  level  of  familiarity  with  system  operation 
was  gained  and  required  changes  in  design  and/or  operational 
procedures  had  been  determined,  a  more-permanent,  all-weather 
installation  could  be  planned. 

It  is  anticipated  that  the  optical  head  of  the  LIDAR 
system  could  be  mounted  on  an  existing  or  easily  installed 
deck-mounted  pedestal  such  as  used  to  support  shipboard 
signaling  lights.  Utilization  of  an  existing,  manually 
aimed  pedestal  should  simplify  shipboard  installation  and 
should  be  accompanied  by  some  reduction  in  cost. 

Operation  of  the  shipboard  LIDAR  could  be  eventually 
performed  by  one  person,  probably  the  ship’s  meteorological 
officer.  A  team  of  one  or  two  scientists  working  in  con¬ 
junction  with  the  ship's  officer  would  be  required  to 
initially  operate,  test  and  evaluate  the  system's  performance 
and  suitability  for  shipboard  use.  Eventually  the  system 
would  be  operated  by  the  ship's  meteorological  officer  and 
used  to  collect  data  on  operational  cruises,  thereby 
generating  the  necessary  data  base  of  marine  aerosol  ex¬ 
tinction  collected  simultaneously  in  the  infrared  wavelength 
bands  of  interest  for  E-0  systems  operation. 

An  optimally  useable  format  for  the  LIDAR  data  should 
be  determined  so  that  it  may  be  readily  used  in  a  PREOS- 
type  FLIR  performance  prediction.  Some  provision  should  be 
made  during  the  shipboard  LIDAR  testing  and  evaluation 
program  to  operate  the  system  along  with  a  FLIR  performance 
prediction  scheme  such  as  PREOS  and  simultaneously  with 
actual  FLIR  measurements  performed  under  well  documented  test 
conditions.  In  this  way  the  correspondence  of  the  LIDAR- 
PREOS  measurement/prediction  approach  may  be  directly 
evaluated  against  the  performance  of  a  FLIR  system. 
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5 

RECOMMENDATIONS 


In  light  of  the  large  indeterminacy  associated  with 
existing  aerosol  models,  an  important  requirement  exists 
to  develop  a  reliable  shipboard  LIDAR  system  to  be  used 
to  collect  a  near-,  mid-  and  long-wavelength  infrared 
marine  aerosol  extinction  data  base  to  be  used  for  model 
verification,  and  to  provide  a  real-time  measure  of 
infrared  aerosol  extinction  for  use  in  E-0  systems  per¬ 
formance  assessments  at  sea.  Accordingly,  it  is  recommended 
that  a  twofold  program  be  undertaken  to  design,  develop 
and  operate  a  dual -wave length  shipboard  LIDAR  system 
operating  in  the  "eye-safe"  near-  and  far  infrared  wave¬ 
length  bands.  This  system  would  be  used  on  several  cruises 
to  collect  a  marine  aerosol  extinction  data  base  and  to 
develop  a  prototype  instrument  for  widespread  shipboard 
deployment.  The  latter  will  be  used  to  provide  routine, 
in-situ  measurements  of  marine  aerosol  effects  on  Navy 
E-O/IR  systems. 

The  specific  approach  recommended  in  carrying  out  the 
shipboard  LIDAR  development  program  involves  the  design 
and  development  of  an  eye  safe,  1.54  ym  shipboard  LIDAR 
utilizing  proven  MIL-SPEC  modules  together  with  a  sti¬ 
mulated  Raman  scattering  module  to  produce  1.54  ym 
radiation. 

Development  of  the  1.54  ym  source  would  be  based  upon 
a  prior  evaluation  of  candidate  MIL-SPEC  1.06  ym  sources. 
The  two  most  promising  candidates  identified  at  this  time 
are  the  AN/TVQ-2  GLLD  system  and  the  A-6  TRAM  system  as 
discussed  in  section  4.2.  These  units  (presumably  en¬ 
gineering  demonstration  units)  would  be  evaluated  regarding 
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their  availability  for  the  shipboard  LIDAR  program. 
Acquisition  cost  and  time  constraints  would  be  part  of  the 
evaluation.  A  stimulated  Raman  scattering  (SRS)  cell  would 
be  developed  and  interfaced  with  the  1.06  uni  MIL-SPEC  source 
to  produce  an  eye  safe  1.54  ym  operating  wavelength. 

A  second  phase  of  instrumentation  development  is  con¬ 
cerned  with  producing  a  reliable  and  suitably  energetic 
source  of  10.6  ym  radiation.  An  evaluation  of  the  Marconi/ 
Ferranti  system  described  in  section  4.2  together  with 
ruggedized  and  environmentally  conditioned  hybrid-TEA  C02 
lasers  using  an  oscillator-amplif ier  configuration  should 
be  performed.  The  potential  advantages  of  heterodyne 
detection  for  operation  at  C02  wavelengths  should  be  explored 
as  previously  discussed  in  section  4.3.5. 

At  the  same  time  that  first-generation  LIDAR  system 
prototype  hardware  development  is  being  performed,  system 
performance  estimates  and  measurement  planning  to  evaluate 
and  refine  marine  aerosol  models  should  be  undertaken. 
Implementation  of  such  a  plan  would,  in  general  terms, 
involve  the  collection,  reduction  and  analysis  of  at-sea, 
marine  aerosol  extinction  data  gathered  during  several 
cruises  spanning  at  least  a  two-year  period. 

A  parallel  evaluation  and  preliminary  design  effort 
concerned  with  possible  differential-absorption  LIDAR  (DIAL) 
and  laser-Doppler-velocimeter  (LDV)  variations  of  the  basic 
shipboard  LIDAR  should  be  undertaken.  Using  these  additional 
alternative  features,  the  shipboard  LIDAR  could  be  used 
for  humidity  and  wind  profile  measurements  as  well. 

Based  on  the  results  and  experience  gained  during 
the  first-generation,  prototype  system  development  and  de¬ 
ployment,  during  which  the  aerosol-model-validation  data  base 
is  collected,  the  design,  construction  and  testing  of  a 
second-generation,  ruggedized  LIDAR  system  for  eventual 
widespread  shipboard  deployment  should  be  undertaken. 
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DEFINITION  OF  A  REMOTE  SENSING 
SYSTEM  FOR  THE  MEASUREMENT 
OF  MARINE  AEROSOL  EXTINCTION: 
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